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I. INTRODUCTION 

It is well known that three dimensional gravity is an integrable system carrying only 
a finite number of degrees of freedom [1, 2]. Despite its apparent simplicity, this theory 
has not been yet fully solved, especially at the quantum level. Such a solution would be 
relevant to people interested in quantum gravity since 3 dimensional quantum gravity raises 
many of the key issues involved in the quantization of gravity such as: the problem of time, 
the problem of the construction of physically relevant observables, the coupling of quantum 
gravity to particles and matter fields, the emergence of semi-classical space-time geometry, 
the fate of black holes in a quantum geometry, the effect of cosmological constant, the 
holography principle, the quantum causal structure, the role of diffeomorphisms, the sum 
over topologies... The notable exception is the absence of gravitons in three dimensions. 

The first problem one faces when dealing with the quantum theory is the variety of tech- 
niques and strategies one can use in order to quantize the theory [3]. The most prominent 
ones being the ADM quantization, the Chern-Simons quantization and the spin foam ap- 
proach. ADM formulation has the virtue to be expressed in geometrical variables but this 
quantization scheme has not been very successful except in a few simple cases [4]. The 
Chern-Simons quantization is the most successful quantization scheme so far, which has 
allowed the inclusion of a cosmological constant and particles [5-7]. The drawback of this 
approach is that there is a clear lack of geometrical and physical understanding when work- 
ing with the Chern-Simons observables. Moreover, this quantization scheme is tailored to 
three dimensions and the techniques developed there have no chances to be exported to the 
higher dimensional case. Eventually, the spin foam approach is the most promising approach 
in a direct quantization of higher dimensional gravity [8, 9]. Being a space time version of 
loop quantum gravity it has an interpretation in terms of quantum geometry [10]. Indeed, 
the first model of quantum gravity to be ever proposed by Ponzano and Regge in 1968 was 
a spin foam model for Euclidean quantum gravity without cosmological constant [11, 12]. 

Among several strategies, coupling gravity to particles is the simplest and most effective 
way to understand the physics and dynamics of quantum spacetime and the construction 
of observables. At the classical level, there is an extensive literature on this subject using 
a diversity of related descriptions: the polygonal t'Hooft approach [13, 14], the Waelbroeck 
discrete phase space approach [15-17], the conformal gauge approach [18], the Chern-Simons 
formalism [19, 20] and the continuum first order formalism [21, 22]. Almost all these works 
concentrate on the case of spinless particles. At the quantum level, most of the explicit results 
concerning the scattering of particles have been obtained using semi-classical techniques 
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[23, 24] or the Chern- Simons formulation [25]. 

In our paper, we consider the spin foam quantization of three dimensional gravity coupled 
to quantum interacting spinning particles. We revisit the original Ponzano-Regge model in 
the light of recent developments and we propose the first key steps toward a full understand- 
ing of 3d quantum gravity in this context, especially concerning the issue of symmetries and 
the inclusion of interacting spinning particles. The first motivation is to propose a quanti- 
zation scheme and develop techniques that could be exported to the quantization of higher 
dimensional gravity. As we will see, the inclusion of spinless particles is remarkably simple 
and natural in this context and allows us to compute quantum scattering amplitudes. This 
approach goes far beyond what was previously done in this context by allowing us to deal 
with the interaction of particles. The inclusion of spinning particles is also achieved. The 
structure is more complicated but the operators needed to introduce spinning particles show 
a clear and beautiful link with the theory of Feynman diagrams [26]. 

Another motivation of our work is to give a unified picture of possible quantization 
schemes. First, we show that the spin foam quantization is related to the discrete Waelbroeck 
hamiltonian approach which we generalize to include spinning particles. This approach is 
known to be related to 't Hooft polygonal approach [17] and this gives a first step towards the 
understanding of scattering amplitude in the 't Hooft approach as a spin foam model. We are 
aware of an independent work to come [27] in the context of loop quantum gravity leading to 
results similar to ours. We also completely unravel in our work, and in a companion paper 
[28] the link of the Ponzano-Regge model with the Chern-Simons quantization. We give a 
complete treatment of the gauge symmetries of the Ponzano-Regge model generalizing the 
work done in [29]. This opens the way to a finite quantization of Lorentzian gravity in the 
spin foam approach [30-32]. We introduce the notion of Ponzano-Regge observables and 
show that the computation of the expectation values of these PR-observables leads to gauge 
fixing, introduction of time or inclusion of particles, depending on the choice. 

We will proceed as follows: In section (II), we review three dimensional gravity in the 
first order formalism, the description of spinning particles in this context and some facts 
about path integral quantization and transition amplitudes. In section (III), we review 
the Ponzano-Regge model, its gauge symmetries and their gauge fixing. In section(IV), 
we describe the notion of Ponzano-Regge observables, their application to gauge fixing, 
or insertion of particles and some duality properties as well as the link of Ponzano-Regge 
quantization with Chern-Simons quantization. In section (V) we present a hamiltonian 
description of three dimensional gravity allowing the inclusion of spinning particles and its 
quantization. In section (VI), we present the insertion of interacting spinless particles. In 
section(VII), we describe the inclusion of interacting spinning particles and the notion of 
particle graph functional. In section (VIII), we compute explicitly, starting from our general 
definition the physical scalar product and the action of he braid group. We show that 
we recover what is expected from canonical quantization and that the braid group action is 
governed by the quantum group D(SU (2)), known as the kappa deformation of the Poincare 
group. 
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II. 3D CLASSICAL AND QUANTUM GRAVITY AND PARTICLE COUPLING 



A. 3D classical gravity 

In this part we recall briefly the properties of 3dimensional Euclidean gravity in the first 
order formalism [9]. 

We consider the first order formalism for 3d gravity. The field variables are the triad 
frame field (i — 1, 2, 3) and the spin connection uj 1 ^. The metric is reconstructed as usual 
from the triad g^ v = e l rjijel where r\ = (+, +, +) for Euclidean gravity. In the following, we 
will denote by e\uj l the one-forms e^dx^ , uj^dx^ . We also introduce the SU(2) Lie algebra 
generator Jj, taken to be i times the Pauli matrices, satisfying [Jj, Jj] = —le^k f] kl Ji, where 
eijk is the antisymmetric tensor. The trace is such that tr(JjJ,) = —2Sij. One defines the 
Lie algebra valued one-forms e = e* Jj and uj = uj l Ji. The action is 

where A is the antisymmetric product of forms and F{uj) = duj + uj A uj is the curvature of 
uj. The equations of motion of this theory are 

d u e = 0, F(w) = 0, (2) 

where d^ — d + [uj, ■} denotes the covariant derivative. If M possess some boundaries dM, 
the variation of the action is not zero on on-shell configurations but contains a boundary 
contribution 

SS =^qI tr(eAM- (3) 



dM 

This boundary term vanishes if one fixes the value of the connection on dM. The gauge 
symmetries of the continuum action (1) are the Lorentz gauge symmetry 

uj -> g~ 1 dg + g~ l ug, (4) 
e -> g^eg, (5) 

locally parameterized by a group element g, and the translational symmetry locally param- 
eterized by a Lie algebra element <fi 

uj — > uj, (6) 
e -> e + <1,0. (7) 

and which holds due to the Bianchi. identity d^F = 0. This supposes that <fi = on 
dM. The infinitesimal diffeomorphism symmetry is equivalent on-shell to these symmetry 
when we restrict to non-degenerate configurations det(e) ^ 0. The action of an infinitesimal 
diffeomorphism generated by a vector field £ M can be expressed as the combination of an 
infinitesimal Lorentz transformation with parameter uj^ and a translational symmetry 
with parameter e M £ M . The conjugate phase space variables are the pull-back of (uj, e) on a 
two dimensional spacelike surface, their Poisson brackets being 

K,ej} = ^V (8) 

The generator of the translational gauge symmetry is given by the pull-back of the curvature 
on the two dimensional slice, whereas the pull-back of the torsion generates the Lorentz gauge 
symmetry. 
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B. Classical particles 



We will be interested in the coupling of gravity to particles. For the reader interested in 
a more detailed and complete treatment of this problem at the classical level we recommend 
the reference [24] which contains a comprehensive analysis of this problem. It is well know 
[2] that the metric associated with a spinning particle of mass m and spin s coupled to 3 
dimensional Euclidean gravity is a spinning cone 

ds 2 = (dt + AGsd^f + dr 2 + (1 - AGmfr 2 dtp 2 , (9) 

where m is the mass of the particle and s is its spin 1 . This is a locally flat space, t is the 
Euclidean time coordinate, r the radial coordinate measuring at fixed time, the geodesic 
distance from the location of the particle along constant tp geodesies and tp is an angular 
coordinate with the identification tp — > tp + 2n. When 4Gm < 1 this space can be identified 
with a portion of Minkowski space. Lets consider the wedge < tp < 2tt(1 — AGm), the 
spinning cone is obtained after an identification of the two faces of the wedge by a translation 
along the t axis of length 8ttGs. Around r = 0, which is the location of the particle, there 
is a deficit angle of 8irGm and a time offset of length 87iGs. The mass of the particle is 
necessarily bounded by 1/4(7. A frame field for this geometry can be given by 

e = J dt + (cos</?Ji + sin tpj 2 )dr + ((1 — AGm)r(costpJ 2 — siiK^Jx) + AGsJ ) dtp, (10) 

and the spin connection by 

tij = 2GmJ dtp. (11) 

The torsion and the curvature have a distributional contribution at the location of the 
particle 2 

d w e = 8irGsJ 5 2 (x)d 2 x, (12) 

F(u) = A7iGmJ 5 2 (x)d 2 x, (13) 

where the delta function is along the plane t = cste. Since the torsion is the generator of 
Lorentz gauge symmetry we see that having a spin means that this symmetry is broken at 
the location of the particle, also the mass is breaking the translational symmetry at the 
location of the particle. We can explicitly see that this is the case if we perform a Lorentz 
transformation labelled by g^ 1 and then a translational transformation labelled by —0, the 
equations (12) then become 



d w e = ATvGj5 2 (x)d 2 x, (14) 
F(u) = 47iGp5 2 (x)d 2 x. (15) 



where j, p are the following Lie algebra elements 



p = mgj g \ (16) 
j = 2sgJ Q g~ 1 -m[gj g'\(f)]. (17) 



1 in H unit so it has dimension of an angular momenta. 

2 we use the distributional identity ddip = 2n5 2 (x)dxdy. 
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p is the momenta of the particle and j the total angular momentum, they satisfy the con- 
straints 

— -tip 2 = m 2 ; — -tr(pj) = 2ms. (18) 

From the canonical point of view these constraints are first class [24], the mass constraint 
generates time reparameterization and the spin constraints generates U(l) gauge transfor- 
mation g — > gh. Due to the breaking of symmetry at the location of the particle the gauge 
degrees of freedom g, <p become dynamical (modulo the remnant reparameterization plus 
U(l) gauge symmetry): g describes the Lorentz frame of the particle, 4> describes the po- 
sition of the particle. Moreover the knowledge of p, j is enough to reconstruct g,4> modulo 
the remnant gauge symmetries. Indeed gH is determined by (16). If we denote by x a ± the 
position of the particle perpendicular to the momenta, = '^jrP + x±, then x± = ^[j,£>], 
also j 2 = s 2 + m 2 x\. 

We can easily understand the canonical commutation relations of p, j from the equations 
(14,15). Since the LHS of (14) is the generator of Lorentz transformations and the LHS of 
(15) is the generator of translational symmetry, these constraints are first class and from 
their canonical algebra we can easily deduce that the Poisson algebra of p, j is given by 

{f, /} = -2e ab %, {f, p b } = -2e abc p c , {p a , p b } = 0. (19) 

This analysis shows that, instead of treating the gravity degrees of freedom and the 
particle degrees of freedom as separate entities, we can reverse the logic and consider that 
the equations (12,13) are defining equations for a spinning particle. This allows us to describe 
a particle as a singular configuration of the gravitational field giving a realization of matter 
from geometry. The 'would-be gauge' degree of freedom [33] are promoted to dynamical 
degrees of freedom at the location of the particle. This is the point of view we are going to 
take in this paper. In order to get the equations (12,13) from an action principle we have 
to add to the gravity action (1) the following terms 

S Pmt3 (e,uj) = ~\j ' dt tr[(me t + 2su t )J ], (20) 

where the integral is along the worldline of the particle. This action describes a 'frozen' 
particle without degrees of freedom. We have seen in (14,15) that the particle degrees of 
freedom are encoded in the former gauge degrees of freedom. To incorporate the dynamics 
of the particle we perform the transformation 

w -> to = g _1 ug + g^dg, e -> e = g^(e + d w <p)g (21) 

the action (20) becomes 

Sp m ,Ae,u) = ~ J dttY{e t p + u t j] + S Pmia (g,<P) (22) 

where the first term describes the interaction between the particle degree of freedom and 
gravity. The second term 

Sp m , 3 (9,<f>) = ~\m J dtti(g~ l <PgJ ) - s J dtti(g~ l gJ ), (23) 
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is the action for a relativistic spinning particle in a form first describe by Sousa Gerbert [24] . 
One sees again that the original gauge degree of freedom are now promoted to dynamical 
degree of freedom describing the propagation of a particle. 
Note that S Pm _ s (g,(f)) = S Prns (ge, -0) where 




so this action describes a particle carrying both spin s and — s, which is necessary in order 
to have P, T invar iance. 



C. Quantum amplitudes 

So far we have described the classical features of 3d Euclidean gravity in the first order 
formalism. We want now to present some general features of partition functions, transition 
amplitudes and particle insertions which are necessary in order to deal with the explicit 
quantification of the theory that we will present in the next sections. Before going on, one 
should keep in mind that the theory we are studying is classically equivalent to usual gravity 
as long as we restrict to non degenerate configurations of the frame field, det(e) > 0. This 
condition of non degeneracy is not easy to implement at the quantum level and we will 
pursue the quantification without worrying about this issue. We have however to keep in 
mind that the resulting quantum theory is expected to depend on this issue as shown by 
the study of symmetry [34] or by the study of spacetime volume expectation value [35]. 

Let us also emphasize that we are dealing in this paper with the quantification of Eu- 
clidean gravity. Euclidean gravity is a theory which admits an hamiltonian formulation 
and its quantification is a well posed problem as shown for instance by the Chern-Simons 
formulation. Since this is often a source of confusion let us stress that the quantification of 
Euclidean gravity has nothing to do with a Wick rotation of Lorentzian gravity, which has 
never been proven to make sense. This is why in the following we are dealing with quantum 
amplitudes and not ill-defined, Wick-rotated, statistical amplitudes. 

Let us eventually emphasize that the discussion in this section is a formal discussion 
which will help us to introduce all the relevant notions and notations that we will be properly 
defined in the subsequent sections. 



1. Partition function 



At the quantum level, the prime object of interest is the partition function, 
closed manifold M we consider 



J M 



I 



VuVe exp 



IQtiG 



hi 



tr(eAF(w)) 



Given a 



(24) 



In order to have a proper definition of Zm, we should not overcount the configurations 
which are equivalent by gauge symmetries and one should restrict the integral over gauge 
equivalent classes of fields (e, cu) by gauge fixing. We will call 'kinematical observable' a 
general functional of the field 0(e, uj). One of the main object of interest at the quantum 
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level is its expectation value 



Vu)Ve exp 



16nG 



tr(eAF(w)) 



M 



O(e,ou). 



(25) 



A 'physical observable' is a gauge invariant (under Lorentz symmetries and diffeomorphisms) 
observable. Of course, the final interest is in physical observables, but the kinematical observ- 
ables will also be of interest to us. For instance, a gauge fixing procedure or the introduction 
of a particle is realized by the insertion of a kinematical observable. 



2. Transition amplitudes for pure gravity 

In the case of a manifold M with boundary we need to specify boundary data in order to 
have a well defined path integral. The unconstrained phase space of 2 + 1 gravity is given by 
conjugate pairs (e, a;), which are the pull-back of (e, a;) on dM. As we have seen in eq.(3) the 
natural choice is to take the boundary connection ui to be fixed in dM. This amounts to a 
choice of polarization where we consider wave functionals to be functionals of the connection 
^(u)). Another natural choice, is to fix e, hence the geometry, on the boundary. In this case 
we have to add a boundary term to the action 



1 



16nG 



tr(eAw). (26) 

dM 



The path integral, with the set of boundary data, defines for every manifold M with a 
boundary, a wave functional 

G M (u) = [ VuVe e lS[uJ > e] . (27) 

We can choose to split the boundary into initial and final boundary components dM = 
where E/ is equipped with the orientation coming from M and Ej is equipped 
with the opposite orientation. For instance when M = E x /, we have Ej = E x {0}, 
E/ = E x {1}. In this case cu = (u>/,<Dj) and (?M(w/,Wi) is the kernel of the propagator 
allowing us to compute transition amplitudes between two states 



J VuVee^^iOf)®^) = J Vuj <f>* f (u f )G M (u f , ^)^(^). (28) 



This scalar product is interpreted to be the physical scalar product, hence it should be 
positive. It is not, however, a definite positive product: it is expected to have a kernel 
characterized by the hamiltonian constraint. 

A natural basis for gauge invariant functionals of cu is given by spin network functionals 
[36, 37]. Such functionals are constructed from a closed trivalent graph T whose edges e 
are colored by representation j s of SU(2). We denote such states by $(rj g ) (<!>). Such states 
are constructed by first taking the holonomy of the connection along the oriented edges e 
geip) = exp 1 (J.w), then taking the matrix elements of this object Dj-(g E (uj)) in the spin 
je representation and finally by contracting ® g Dj-(g s (u))) G ®eVj E ® Vj s using invariant and 
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normalized tensors C : V± <g> V 2 <g) V3 — > C at the vertices t> of T 3 . So 

$(r,; e -)(^) = (®«C7 S | ® e - D je (g g (u))). (29) 

The construction needs T to be an oriented graph, however the results are independent of 
this orientation in the case of SU(2). Given a graph T, we can consider the vector space 
generated by all spin networks with support T, i-e the space of all cylindrical functions of the 
connection having the form $ = Cj g $(rj g ) where the sum is a finite sum over admissible 
coloring. We denote this vector space by Hr- This vector space can be promoted to an 
Hilbert space with the norm 

\m\ 2 = [ n^e-|<J>G?e-)| 2 , (30) 
J e 

dg being the normalized SU(2) Haar measure. This Hilbert space admits a basis independent 
description 

H r = L 2 (G®/GW), (31) 

where |e| denotes the number of edges of T, \v\ the number of vertices and the action of 
on G' e ' is given by the structure of the graph: g s — > h~^g s h s ( S ), t(e) and s(e) being the 
target and the source of the edge e. 



3. Transition amplitudes with massive spinning particles 

In this paper, we also want to describe the case where particles are present. In this case, 
as we have seen in the previous section, the gauge symmetries are broken at the location 
of the particle, moreover the previously gauge degrees of freedom are now promoted to 
dynamical degrees of freedom of the particle. In order to describe the insertion of particles, 
we have to add the particle action (20) to the gravity action. This breaks the initial gauge 
symmetry. The propagator is given by 

G M {Q)= f VujVee iS[ "> e]+ ZnSpn(«,e). (32) 
Juj\dM=ai 

where n = (m n , s n ) denotes a mass and spin. We can decompose the path integral into an 
integral over the gravitational fields and an integral over the gauge degrees of freedom at 
the location of the particles 

G M ,p n (w,g nii ,g n j)= VcuVe e lS[uJ ' e] Y[G P (ujp n ,ep n ,g n4 ,g n j), (33) 

Jcj\dM=uj n 

G P m A u P^ e P^9i,9f) = I VgV<p exp( I dt - hi[e t p + u t j] + S Pmt ,(g,<j>)). (34) 

Jg\dM=(gi,g f ) J 1 

(up,ep) denotes the value of (ou,e) at the location of the particle. In the integral (33) 
over (u,e), the action is gauge invariant, contrary to the one in (32), and the integration 

3 the space of trivalent intertwiners is one dimensional in the case of SU(2), since we consider only trivalent 
graph we do not need to specify an intertwiner label at vertices. 



9 



is understood to be over gauge equivalence classes. If we concentrate first on the particle 
integral Gp(u,e) and suppose for the sake of the argument that uo = e = along the 
particle worldline, then the interaction term disappears and we are left with Sp ms (g,(p). 
The propagator Gp is then the propagator of a relativistic particle and as such an operator 
acting on the Hilbert space obtained by the quantification of (19,18). Such a quantification 
can be easily described as follows. We take our kinematical Hilbert space to be L 2 (G), 
such an Hilbert space is spanned by the Wigner functions D^ k (g) where the spin / is an 
half integer and n, k are the representation indices. We define the operators p to act by 
multiplication, and j to act as the right invariant derivative 

p a DUg) = mtiigJog-^D^ig); jaD^g) = -D^g). (35) 

The mass constraint p 2 = m 2 is trivially satisfied, the spin constraint implies that k = — s 

p^jaD^g) = -mtiigJog-iJ^D^Jag) = mD^gJo) = 2msD I n _ s (g). (36) 

The physical Hilbert space associated with the particle is the usual Poincare representation, 

n m , s = V I = {D I n _ s (g)\I-seN, \n\<I}. (37) 
i\i-seN 

Instead of labelling the particle propagator by couples g^gf we can label it by a pair of 
Lorentz indices h,If- This means that the particle propagator Gp ms j i j f (ou,e) is viewed 
as an operator in Hom(V r ., Vj / ). The restriction we made on the value of (u>,e) along the 
particle world-line can be relaxed, it doesn't change our argument but just the value of 
Gp m 3i 7 i: / r This shows, in the case of particle insertions, that the propagator depends not 
only on a given interpolating manifold M but also on additional data characterizing the 
evolution of the particles. The labels Ii,If are labels of the total angular momenta, as 
we have seen, we can also understand them as labelling the position of the particle in the 
direction transverse to the momenta since + 1) — + 1) = m 2 [(x 2 L )f — 

We will consider the general case of interacting particles of arbitrary spin. In this case 
the data we need are encoded into what we will call a decorated particle graph and will 
denote by "Id 4 - "I is graph embedded in M such that its open ends v are all lying in DM . 
D = (m e , S e , I Se , It e ,iv) is a decoration of 1 where each edge e of 1 is labelled by a mass and 
a spin (m e , S e ); each starting point of an edge e is denoted by s e and labelled by a SU(2) 
representation V} , each terminal point of an edge e is denoted by t e and labelled by a SU(2) 
representation Vi te ; finally each internal vertex v is labelled by an intertwiner ^ contracting 
the Lorentz representations I (see figure 1). 

Given a decorated manifold with boundary M,~[ D we now want to describe the construc- 
tion of the kinematical spin network states and their transition amplitudes. We have to 
consider states that are not gauge invariant at the location of the particle. Such kinematical 
states are described by spin networks with open ends [38]. We denote by v the internal 
vertices of T and by v the open ends of T, they are identified with the open ends of the par- 
ticle graph 1. We also denote by e the internal edges of T and by e v the edges of T ending 



4 We could equivalcntly call these graphs Feynman graphs since they label all possible Feynman graphs of 
3d field theories, however this terminology usually suppose that there is a specific field theory behind the 
construction of this graphs and that we should sum over them. It is an interesting problem to describe 
such a theory but this is not the goal of this paper 
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FIG. 1: boundary and particle graphs 



on an open end v. Given a graph T whose internal edges are colored by j g and external 
edges colored by I v , we can construct in the same spirit as before a spin network functional 
*&(T,jg,i v )(&)- The difference with the previous case is that now $(r ,j g ,i v ) is not a scalar 
but take value in <8> v Vi v . The construction is similar to the previous one, let us denote by 
ge{oj) the holonomy along internal edges and g v {ui) the holonomy along external edges. The 
spin network functional is given by 

®(T, je ,i v )(u) = ®e D je (g g (Q)) ® v D Iv { 9v {uj))) G ® v V Iv . (38) 

The intertwiners C„ live only at internal vertices where a pairing occur, there is no con- 
traction at external vertices. As before, we can consider the linear combination of all spin 
networks supported on a graph and define the basis independent spin network Hilbert space 

Hr,i v = L 2 (((jU+M -> ® v V Iv )/GM) , (39) 

where |e| denotes the number of internal edges of r,|i>| the number of open ends of T, \v\ 
the number of internal vertices and the action of on G^ is given by the structure of the 
graph. This space is our kinematical Hilbert space. 

It should be clear, from the discussion of the classical particle, that the spin label I v on 
open edges have the meaning of the total angular momenta carried by the particle. The 
individual vectors in the representation Vj v label the direction of the angular momenta. I v 
doesn't have the meaning of the spin of the particle, which is the component of the angular 
momenta along the direction of the particle momenta, this is important to stress since this 
is often a source of confusion. The spin and the mass appear dynamically in the choice of 
the propagator determining the physical scalar product. 

Due to Lorentz gauge invariance not all choices of I v are admissible and when a choice is 
admissible not all vectors in ® v Vi v appear in the image of $(r,j g ,j„)- The image of ^rr,j g ,i v ) 
lies in the invariant subspace (<S> v Vj v ) G which is the space of intertwiners C — > <2) v Vi v , where 
C denotes the trivial representation of SU(2). 

The full propagator is not a scalar function but an operator depending on the decorated 
particle graph: 

Gm^Q/M) ■ {®vVi v f - (®v f V Ivf ) G . (40) 
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The amplitude is given by 



($/l^>M,n D = J VQ^fiQf^GM^iQ^Qi^iiQi)). (41) 

The measure T>uj is the Ashtekar-Lewandowski diffeomorphism invariant measure [37]. Due 
to diffeomorphism invariance the propagator depends only on the diffeomorphism class of 
the embedding of 1 in M. In the case where M = £ x /, £ being a 2d surface with n 
punctures, the physical scalar product on Hr,i v is given by the following decorated graph: 
1 is the union of unbraided oriented segments e v joining Vi with Vf. The coloring of each 
line by (m ev , S ev ) labels the mass and spin of the particle, also I s {e v ) = ^ and h{e v ) = Iv r 



III. THE PONZANO-REGGE MODEL AND GAUGE FIXING 



The Ponzano-Regge model is constructed from the continuum partition function 



'M 



VuVe exp 



16nG 



tr(eAF(w)) 



M 



(42) 



by considering a triangulation A (and its dual A*) of M, and replacing the set of configura- 
tions variables by discrete analogs in the spirit of lattice gauge theory. The connection field 
is replaced by group elements g e * associated to the dual edges e* of A*, and representing 
the holonomy of the connection field along these edges. The frame field is replaced by Lie 
algebra elements X e associated to the edges e of A and representing the integration of e 
along these edges. The curvature 2-form is now represented as group elements G e living on 
the edges e (or dual faces /*), and obtained as the ordered product of the group elements 
g e * for dual edges e* C /*, upon the choice of a starting dual vertex on the dual face. 
• The discretized partition function is thus 



Zpr(A) = [ J] f dg e *) (n 

V e . ^SU(2) J V e 



dX e 1 exp 



su(2) 



^Ttr(X e G e ) 



(43) 



One can then integrate over the X e variables (see appendix (B)), 



(44) 



where 5(g) is the delta function over the group 5 . The Ponzano-Regge partition function is 
recovered from (44) by expanding the 5 function using the Plancherel decomposition, and 
then integrating over g e * using recoupling identities for SU(2) [8]. The result is a realization 
of (24) and can be expressed as a summation over coloring of a product of 6j symbols. 



Z PR (A) 



{je} e t 



Jet 1 Jet 2 Jet 3 
Jet 4 Jet 5 Jet 6 



(45) 



5 As it is explained in the appendix (B) it is the delta function on SO (3) and not SU(2). It is also explained 
how we can modify (43) in order to get the delta function on SU(2). 
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where the summation is over all edges of A and the product of 6j symbols is over all 
tetrahedra. dj = 2j + 1 denotes the dimension of the spin j representation and denotes 
the six edges belonging to the tetrahedra t. j e plays the role of a duiscrete X e , it is therefore 
interpreted as a length and the Ponzano-Regge sum is a sum over all geometries supported 
by a triangulation. 

• In the case of manifold with boundaries but no particles, we can define the physical 
scalar product between spin network states ®r f ,j Bf e 7~lr f , ^Ti,j Si £ 7^r\ hi a similar way 

< 'lY,.,.. 'IY..,.. >pr= [ J] / d 9e* ) n S ( G e)®r f , jgf (//... H'r..,. (9*), (46) 

the integral being over all group elements associated with dual edges including the ones e/, 
lying in the boundary. As before we can recover Ponzano-Regge like expressions giving the 
physical scalar product (28) as a state sum model 

< >™= sii4.nl % % % } , (47) 

{j e } e t < J * 4 J *5 J6t 6 J 

where the summation is over all edges of A which do not belong to the boundary, the coloring 
of the boundary edges e/, is fixed to be j Sf , je v 

• In the general case of manifold with boundaries and with particles we can also define 
a transition amplitude which gives an explicit realisation of (41). This one has never been 
written before 6 and the explicit description of it is one of the main result of this paper. We 
will present its construction in detail in section VI. In the case of interacting spinless particles 
the result is very simple. The particle graph 1 is living on the edges of the triangulation A, 
and we choose a decoration D where all spins and angular momenta are 0. In this case the 
amplitude including boundaries and particles is given by 

<*,m»» - e (nt) (n*^>) n { t: t i: } • < 48 » 

the summation is over all edges of A which do not belong to the boundary, the coloring of 
the boundary edges ef, is fixed to be je f ,jei- The edges in 1 carry a factor Xj(h$), which 
is the trace in the representation of spin j of the group element h e = exp(6'J ) (Al). 9 is 
half the deficit angle created by the presence of a mass m, 9 — AirGm. 

The general case involving interacting spinning particles is similar but more involved 
and described in detail in section VII. For all these amplitudes we can recover Ponzano- 
Regge like expressions as state sum models. However the naive definitions of these models 
is plagued with divergences. We need to do a careful analysis of the gauge symmetries of 
the discrete model in order to understand these divergencies and to take care of them. This 
is what we are presenting now. 



A. The gauge symmetries of the Ponzano-Regge model 

In this part we present the identification of the gauge symmetries of the Ponzano-Regge 
model, and their gauge fixing. We discuss also the case of a manifold with boundaries and 



some hints concerning the insertion of spinless particles were already given in [39] 
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insertion of particles. This work completes the work begun in [29] where we only gauge fixed 
one of the symmetries. Note that we corrected here the exact form of the action (43) and 
as a consequence, the formula are slightly different from [29]. 

1. Discrete symmetries of the Ponzano-Regge model 

The gauge symmetries of the continuum action (Lorentz and translation gauge symme- 
tries) have been described in the first section. The discrete analog of the Lorentz symmetry 
is parameterized by group elements k v * living at the dual vertices of the triangulation. It 

acts as 

9e* -> K e ]g e *k Sef (49) 
G e - K t ]G e k ste (50) 
X e - K t \X e k ste (51) 

where s e * and t e * denote the dual vertices source and target of e*, st e denotes the dual vertex 
which is the starting point for computing the curvature on the dual face /* ~ e. 

To express the discrete analog of the translational symmetry, we observe that the discrete 
action appearing in (43) can be written as 

S = tr(^X e P e ) (52) 

e 

where P e = P^J a is the projection of G e on the Lie algebra and is defined by 

G e = ujd + P a e J a , (53) 

where u 2 e + P^P ea — 1. In the continuum we have seen that the translational symmetry is 
due to the existence of Bianchi identity. At the discrete level we also have a Bianchi identity 
for each vertex of the triangulation: [40] 

XliKr'Gf^K = Id, (54) 

eDv 

where the product is over all edges meeting at v, k e v are group elements dependent on g e * 
characterizing the parallel transport from a fixed vertex to st e and e(e,v) = ±1 depending 
on whether e is pointing toward v or in the opposite direction. We can write this identity 
in terms of the P e . One finds that there exist Lie algebra elements fig and scalars U v e which 
can be expressed in terms of the P e 's for edges meeting at v, such that 

^TeM)^)- 1 (u v e p e + [ni,p e ]) k = o. (55) 

It is easy to check that the discrete analog of the translational symmetry is then parame- 
terized by Lie algebra elements Q v living at the vertices of the triangulation and acts as 

X e - X e + e(e, v) {U: {K^K)' 1 ) - [^, (A^(^)- 1 )]) (56) 

for v C e. 
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2. Gauge fixing of the symmetries 

Having identify the discrete gauge symmetries of the discrete action (43), it has been 
observed [29] that the infinite gauge volume of the translational symmetry was actually 
responsible for the divergencies of the Ponzano-Regge model. We then proposed a procedure 
to gauge fix this symmetry. The Lorentz gauge group being compact in euclidian case, we 
ignored the gauge fixing of the Lorentz symmetry and concentrated on the translational 
symmetry. Having in mind the possible applications for the Lorentzian case where the 
Lorentz group is no longer compact, we now need to be able to perform both gauge fixing 
at the same time in a consistent way. This is the purpose of this part. 

Both discrete gauge symmetries act at the vertices of a graph (1-skeleton and dual 1- 
skeleton of A), on variables living on the edges of this graph. The usual method to gauge 
fix such a symmetry is to choose a maximal tree, i.e. a connected subgraph touching every 
vertex without forming a loop, then to use the gauge symmetries at the vertices to gauge fix 
all the variables living on the edges of the tree. This method has been applied to the Lorentz 
symmetry as a gauge fixing for the non-compact spin-networks [41], and to the translational 
symmetry for the Ponzano-Regge model [29]. Recall that this procedure uses the action of 
the gauge symmetries at every vertex except one, taken as the root of the tree. This gives 
rise to a remaining global symmetry that has to be studied in a second time. 

• We want to perform both gauge fixing. We thus choose a maximal tree T in the 1- 
skeleton of the triangulation A and a maximal tree T* in its dual 1-skeleton. If \v\ and \t\ are 
the number of vertices in the 1-skeleton and dual 1-skeleton, these trees contain respectively 
\v | — 1 and \t\ — 1 edges. We use the gauge symmetries to fix 

g e * = 1, Ve* G T*, (57) 
X e = 0, Ve G T. (58) 

We now need to identify the corresponding Fadeev-Popov determinant. The jacobian for 
fixing the translational symmetry has been computed in [29]. Using the fact that T* is a 
tree, one can show that the jacobian associated to the gauge fixing of the Lorentz symmetry 
is I. Finally we have for the Fadeev-Popov determinant 

Ai,p = n(( f/ e e ) 2 +l^e e | 2 )|^ e |- (59) 

eeT 

Using the argument given in [29] we can prove that this determinant is actually 1 while 
included in the partition function. This is due to the fact that the product of S(G e ) for 
e ^ T is enough to impose P e = on every edge as long as T does not touch the boundary 
or the particle. The net result is that there is no Fadeev-Popov appearing in the partition 
function after the double gauge fixing. After this gauge fixing procedure we have 

Z[A\{T,T*}] = ( II / d sA H*(G e ). (60) 

\e*^T* ^su(2) J e ^ T 

• After this gauge fixing procedure, we are left with a triangulation which possesses only 
one vertex and one dual vertex. However, there are still some invariances associated. The 
remaining Lorentz invariance acting at the unique dual vertex is a diagonal AdG invariance 
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g e * — ■> ^g./c. It has been shown in [41] that this remaining invariance can be gauge fixed 
using the measure 

dfi(gi, g N ) = d/j,(g 1: g 2 )dg 3 ...dg N , (61) 
where dg is the Haar measure on G and d/i(gi,g 2 ) is defined by 

/ d/i(g 1 ,g 2 )f(g 1 ,g 2 ) = dhdx, f(h,s(x)) (62) 

Jg 2 JHxG/H 

where H is the Cartan subgroup of SU(2) and s : G/H — > G is a given section 7 . This choice 
of measure allows to gauge fix the remaining Lorentz invariance. 

Let us now consider the case of the remaining translational symmetry. This symmetry is 
supposed to arise from the Bianchi identity around the last vertex. As this vertex is the last 
one, every edge starts and ends on it. The corresponding Bianchi identity (54) thus involves 
both G e and G?" 1 . The key point is to understand the order of the elements in (54). This 
order is related to the topology of the neighborhood of the vertex. Since we started from 
a triangulation of a manifold, the neighborhood of each original vertex was a three-sphere. 
Now we remove the edges corresponding to a tree, which by definition has no loops. Hence 
its tubular neighborhood has the topology of a sphere and we are left with a last vertex 
whose neighborhood is still a 3-sphere. The Bianchi identity is trivial which means that 
there is no residual action of the translational symmetry. This just means that our original 
parameterization of the gauge group was redundant, and that the gauge group was actually 
su(2) |,;| -'. 

3. Gauge fixing: a mathematical argument 

There is a beautiful and simple mathematical argument which allows to understand the 
necessity of gauge fixing which goes as follows. We know that in the theory of integration 
what we should integrate over a manifold of dimension n is not a function but a density of 
weight n. A delta function on the group SU(2) is a distribution, that is a density of weight 3, 
i-e the only thing that make sense is the quantity dgS(g) and not 5(g). Before gauge fixing, 
the integration manifold is a manifold of dimension |/|, the number of faces of A and the 
integrand is a density of weight |e|, the number of edges of A. The dimension of the space 
does not match the density weight of the integrant and the integration does not really make 
sense. After gauge fixing the density weight of the integrant (60) is 3(|e| — \v\ + 1) and the 
dimension of the integration manifold is 3(|/| — \t\ + 1), where \v\ (resp. \t\) denotes the 
number of vertices (resp. tetrahedra) of A. The difference between the two is given by 

-3 X = 3(|t|-|/| + |e|-M), 

where x is the Euler characteristic of the 3d manifold M. If M is a closed manifold its Euler 
characteristic is zero, therefore the integral (60) is a priori well defined. If M is not closed 
x(M) = l/2x(dM) = (1 — g), where g is the genus of the boundary if dM is connected. 
Therefore if the boundary does not contain any sphere or torus — 3%(M) = 3(g — l) is strictly 



7 We present here only the measure in the case of a unique Cartan subgroup, relevant for SU(2). The case 
of many Cartan subgroup relevant for SL(2,R) can be treated in the same way, see [41] 
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positive. The result of integration in this case is a density of positive weight 3(g — 1) which 
is exactly half the dimension of the moduli space of flat SU(2) connection on DM. This is 
the correct result since we will see that in the case of manifold with boundaries, the quantity 
(60) should be interpreted as the integration kernel for quantum transition amplitudes. One 
should be careful however that the counting argument we have just given is naive and does 
not always guarantee the defmiteness of the integral (60) 8 . 

In terms of the mathematical argument given at the beginning of this part, the result of 
this integration can be evaluated to be 



A e *g T *dg e 



f\ e ^T<lG e 



(63) 

G e =0 



where A denotes the wedge product of forms. When the manifold is closed this is just a 
jacobian, when the manifold admits boundaries this is a form of positive degree by the 
previous argument. This argument does not guarantee the defmiteness of this jacobian 
since there still could be situations in which (g e *) e *eA* does not depend on some component 
of (<7 e ) ee A- 111 this case the jacobian will be infinite. We will analyze more precisely the 
finiteness issues in section VIII D. 



4- Inclusion of boundaries and particles 

We now discuss the generalization of the gauge-fixing procedure to the case of a manifold 
with boundaries, when we try to compute transition amplitudes between spin-networks 
states. We have seen that at the classical level the Lorentz symmetry also acts on the 
boundary, while the translational symmetry is such that the parameter of transformation 
has to vanish on the boundary. This is translated into the fact that the physical scalar 
product between spin-network states 

< 'IV.,.. 1^ >PR= ( II / ^e* ] II <I>r.,„, (9^r t , hi (<fcj, (64) 

has a Lorentz symmetry acting at all the dual vertices of the triangulation (including the 
dual vertices lying on the boundary) while the translational symmetry is restricted to the 
bulk vertices. The gauge fixing of such amplitudes involves then a maximal tree T of the 
triangulation, touching at most one vertex on the boundary (recall that there is one vertex 
in the maximal tree that we actually don't use for the gauge fixing). The gauge fixing of the 
Lorentz symmetry involves a dual maximal tree T* touching every dual vertex, including 
those on the boundary. 

As discussed in the classical case, it is expected that the introduction of a particle will 
break a part of the gauge symmetries of gravity, turning them into particle degrees of freedom 
and modify the discussion of the gauge fixing. Actually, the inclusion of a particle on an 
edge e breaks the translational gauge symmetry as does the inclusion of mass classically. 
As we will discuss explicitly later, the discrete transformation (56) is not a symmetry of 
the discrete action anymore. A maximal tree T for the gauge fixing of the translational 
symmetry has to be chosen outside the vertices lying on the trajectory of the particle. 



In a similar way as having a negative naive superficial divergences degree of a Feynman integral does not 
always guarantee its convergence. 
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B. BRST invariance 



We have computed the gauge fixed partition function in terms of a gauge fixing procedure 
involving two maximal trees. In a general gauge theory the gauge symmetry is obviously no 
longer visible after gauge fixing. However there is an in-print of this symmetry in the gauge 
fixed theory which is BRST invariance assuring us that the gauge fixed partition function is 
independent of the gauge choice. In our case we have to check if the gauge fixed partition 
function 

Z[A\{T,T*}] = ( J] I dg e A l[5(G e ), (65) 

depends on the choice of (T, T*). It is shown in a companion paper [28] that the gauge fixing 
procedure is independent of the choice of maximal trees i.e. 

Z[A\{T U T*}] = Z[A\{T 2 ,T*}\. (66) 

This property is a remnant of the original gauge + translational symmetry. It follows from 
the invariance properties of the path-integral but is proved more easily using the graphical 
tools introduced in [28]. It is also proved in this paper that the partition function does not 
depend on the choice of the triangulation matching the boundary data so it is a topological 
invariant 

ZlAMTuT*}] = Z[A 2 \{T 2 ,T;}}. (67) 



IV. PONZANO-REGGE OBSERVABLES 

In this section we introduce the notion of Ponzano-Regge observables that are discrete 
analog of 'kinematical observables' presented in section II CI. We show that the gauge 
fixing procedure can be written as the evaluation of the expectation value of a gauge fixing 
operator. We also present the Ponzano-Regge observables realizing the insertion of particles. 
We describe the link between Ponzano-Regge model and Chern-Simons theory and show a 
remarkable duality property between gauge fixing and insertion of particles. 

First, let us remark that if one uses the Plancherel identity 

%)=£|;X,(<?), (68) 

where j is a spin, Xj{9) = ^ T Vj(Dj(g)) is the character of the spin j representation, dj is 
the dimension of the representation and Vq is the volume of the group, then the partition 
function (65) associated with a triangulation A can be written as 

f IP^A(je,Se*), (69) 

where, 

Z PR {A) = Z A ( 3e ,g e *) = lb Xj.(G e ), G e = (70) 
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where e(e, e*) = ±1 is the relative orientation of e* with respect to e. We will call PR- 
observable 9 a function (D(j e ,g e *) and define its expectation value as 

(0)A = J^^l[ e ^ T J dg e *Z A (j e ,g e *)0(j e ,g e *). (71) 

je 

A. Gauge fixing observables 

In particular, an important class of PR-observables are those fixing the lengths along a 
graph T of A 

Or, 3 r(j e ,g e ,) = H S -f^, (72) 

where 8jj> is the Kronecker delta function. And those fixing the holonomy along a dual 
graph T* 

e *er* 

where G [0, 2n] and 80(g) is a delta function fixing g to be in the same conjugacy class as 
h e = exp(26>J ) 10 . More generally, we can define operators fixing edges length along a graph 
T of A and the conjugacy class of g* along a dual graph T* of A*. 

^xr*,^, = n n % w - (74) 

In general the evaluation of the expectation value of ^(r,j r )(r*,e r » ) are Quite involved; however 
if T = T and T* = T* are trees in A and A* then the evaluation simplifies drastically, we 
have the following [28] 

(°(T, 3T){ T*,e i:) )A = (l[d, 7 ) ( n y) Z[A\{T,T*}\, (75) 
Veer / Ve-eT* YH J 

where Vq, Vh denote the volume of G and its Cartan subgroup H . This shows that the 
degrees of freedom living along trees or dual trees decouple from the bulk degrees of freedom. 
We have seen in the precedent section that indeed, such degrees of freedom are pure gauge 
and therefore the expectation value of PR-observables depending only on such degrees of 
freedom is expected to factorize. Equation (75) show that such expectation value does not 
depend for instance on the topology of the manifold M triangulated by A or the decoration 
on the boundary of M. Such PR-observables are therefore not physical observables, they 
couple only to the gauge degrees of freedom and not to the physical degrees of freedom. If we 
consider the observable O r jr where T is not a tree then this PR-observable will couple to the 
dynamical degree of freedom of the theory and it is a physical observable giving us details 
about the geometry of our space-time. Due to gauge invariance two PR-observables O v .r 1 



9 PR stands for Ponzano-Regge of course, these observables are not in general gauge invariant and in this 

case they are 'kinematical observables' 
the normalization and other properties of this function are described in the appendix A 
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O r .r 2 will carry the same physical information if (Ti,^ 1 ) and (r 2 ,j,F 2 ) are diffeomorphic 
colored graphs or if they can be collapsed using a tree to diffeomorphic colored graphs. 

An other interesting evaluation of the expectation value of this operator is when T is the 
one-skeleton of A (in this case we denote it A, by an abuse of notation) and when T* = A* 
is the one skeleton of A*. In this case we introduce a special notation for this object 

ZA(je,0e*) = (C ) (A,j e )(A*,e e ,))A- (76) 

This object does not contain anymore any summation over the length of the edges and 
integration over the conjugacy class of g e *, its evaluation is therefore trivially finite since 

Z A (je, M = / II ^ II XiM e * Ce X e* h Z' e * )x e } ) ( 77 ) 

Jg/h e * eA * eeA 

where dx is the induced measure on G/H and he is defined in (Al). This formula follows 
from the Weyl integration formula (see appendix A). 

B. Chern-Simons formulation 

The remarkable property of this object, which is proven in detail in [28], is that Z&(j e , 8 e *) 
can be interpreted as a quantum group evaluation of a chain mail link La- The link in 
question La is described in detail in [28] but can be sketchily described as follows: One can 
consider the handlebodies H and H* respectively built as the thickening of the 1-skeleton 
and dual 1-skeleton of A. The gluing of H and H* along their boundaries gives a Heegard 
splitting of the manifold M [42]. One can now consider the meridian circles of H and H*, 
and by slightly pushing the meridians of H into H*, one obtains the link La- It should 
be noted that this link is made from two different kinds of components: the ones coming 
from H and the others coming from H*. The quantum group in question here is a kappa 
deformation of the Euclidean Poincare group ISO(3), it is denoted L> K (SU(2)) 11 and can 
constructed as a Drinfeld double of SU(2) [25]. Its relevance in the context of 2+1 gravity 
as been recently analyzed in [43] . The link is colored by spins j for the components coming 
from H and angles 6 for the H* components, j label pure spin representations of the kappa 
Poincare group whereas 6 is interpreted as the mass of a spin zero representation. 

It is well known that, at the classical level, three dimensional gravity can be expressed 
as a Chern-Simons theory where the gauge group is the Poincare group. The Chern-Simons 
connection A can be written in terms of the spin connection uj and the frame field e, A = 
ujiJ % + e,iP l where J 1 are rotation generators and P % translations. Moreover, since the work 
of Witten [44], it is also well known that quantum group evaluation of colored link gives 
a computation of expectation value of Wilson loops in Chern-Simons theory. Our result 
therefore gives an exact relation, at the quantum level, between expectation value in the 
Ponzano-Regge version of three dimensional gravity and the Chern-Simons formulation. 
More precisely, given a link La whose components are colored by j e , 6 e * we can construct a 
Wilson loop functional WL A j e ,e e * {A) by taking the ordered exponential of the Chern-Simons 
connection along the components of La and tracing the result in the Poincare representations 



k = 1 /AG is the Planck mass in three dimensional gravity 
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labelled by j e , 9 e *. If we denote by (-)cs the expectation value in the Poincare Chern-Simons 
theory the result proved in [28] reads 

<<!W)(A*,* e *)}A = z A (j e ,e e *) = (w LA j eA M)) cs . (78) 

Note that the full partition function can be recovered as a sum 

z PR (A) = J2ir[ n^A 2 (M^A(j e A*)- (79) 

U Jh ' W e* 

From the Chern-Simons point of view this gives a surgery presentation of the partition 
function. 



C. Particle observables 

Let us define another type of observables 

&V>.Ue,9s)=l[8s.fi8 9e (G e ), (80) 

where 9 e G [0,27r]. The expectation value of this PR-observable removes (with Sj ej0 ) the 
flatness condition around the edges of 1 imposed by the partition function, and fixes of the 
holonomies around the edges of 1 to be in the conjugacy class of hg e (A5). We will justify 
in greater detail in section VI that the insertion of this operator can be interpreted as the 
presence of a spinless 12 particle of mass m e = Oe/AnG on the edge e. For instance, one 
can prove that the mass is conserved at bivalent vertices of "1. One can also check that the 
expectation value of (80) is given by the expression (48). 

Spinning particles can also be introduced by the insertion of a PR observable described 
precisely in section VII, The structure of the observable is as follows. First we consider a 
decorated particle graph ~~\ D (see section II C 3) living along the edges of the triangulation 
A, we introduce the operator defined for each edge of the particle graph 

6e e ,s e ,i te ,ijje,ge*) = 5 je , [ duAGeUehe^D^iglu^Pf^D 1 ^ ((sl^r 1 ) (81) 

Jg/h 

where 9 e is the deficit angle associated with the mass m e = 9 e /8iiG, S e is the spin, s e ,t e are 
the source and target of the edge e, I Se ,It e are Lorentz indices, D 1 (g) is the representation 
matrix of the group element g in the representation of spin /, Pf t e Is = \h e , S e >< I Se , S e \ 
where \I,S > is the vector of magnetic moment S in the representation I and g% e is the 
group element associated with a path going from st(e) to t* (see section II C 3). Associated 
with each internal vertex of Id we have an intertwiner i% that can be used to contract the 
open indices of Oe e ,s e ,i te ,i Se , to get 

Oln = Yli®viv\ ®e 9e , Se ,I te ,lJ. (82) 
eel 
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more precisely this correspond to a particle which do not carry any angular momenta 
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The expectation of this PR-observable removes (with <5j- e ,o) the flatness condition imposed 
by the partition function around the edges of T, fixes of the holonomies around the edges of 
T to be in the conjugacy class of hg e (A5) and insert a projector forcing the particle to be 
of spin S e . On top of a mass conservation rule there is also a spin conservation rule. In the 
case were all Ps and S"s are zero we recover the operator (80). 

D. Wilson lines and time observable 

Before presenting a detailed justification of the introduction of these massive spinning 
observables corresponding to the inclusion of massive spinning interacting particles it is of 
interest to describe other PR-observables. 

The first one we would like to consider is what we call the identity PR-observable. Let T 
be a graph living in A we denote 

Ir = U 6 ^(G e ). (83) 
eer 

Such an observable removes (with 5j ei0 ) the flatness condition around the edges and insert 
back this condition (with S(G e )). It clearly does not change anything and does not depend 
on T. It is therefore neither physical nor observable, its only interest lies in the fact that 
it appears naturally when we consider the particle observables associated with a zero mass 
and zero spin propagating along T, it is therefore not a surprise that it is trivial. 

More interesting PR-observables are the Wilson observables. Given a graph T* G A* 
whose edges are colored by spins j e * , we can define as in section II C the spin network 
functional $, r *jr*) which depends only on group elements g e *, e* G T*. The Wilson PR- 
observable is given by 

eer* 

Since this observable depends both on the spin connection and the spin on the graph it is 
sensitive to the curvature and metric degrees of freedom. 

The last observable we want to describe here, is by far the most interesting since it allows 
to recover the notion of time at the quantum level. One should however say that since we 
are working in the Euclidean context there is no notion of time orientability and we cannot 
distinguish T and — T, this will be different in the Lorentzian context which will be the 
subject of another work [32]. 

First, we can define the time operator only in the case where the manifold M possesses 
boundaries DM = Ej ]J Ej. In this case we choose L n G A to be a succession of edges of A; 
L = (ei,e 2 , • • -e n ) such that s(e i+1 ) = t(ej) and s(e 1 ) — Vi G Ej, t(e n ) = Vf G Ej. We also 
suppose that L n is not knotted, it is a simple line going from v i to Vf which are two vertices 
of the boundary triangulation. We consider the PR-operator 

n 
i=2 

One can check that the expectation value of this operator does not depend on the number 
of edges composing L n or which edge we choose to be of length j. Ideally we would like to 
define the proper time as the distance between the slices E« and E / along a minimal curve 
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with a fixed initial direction. The insertion of the operator T(j) in the amplitude amounts 
to fix the distance between the two slices but ignoring the information about the direction 
of the curve. It computes a superposition over all initial directions of amplitudes for which 
the distance between the two slices is j. This will be clear in our hamiltonian analysis. 



E. Duality properties 

In this section we would like to present a remarkable duality property between the particle 
PR observables and the gauge fixing PR observables. 

Proposition 1 Consider a graph T in A. The observables O r jr and C*r> e are dual by 
Fourier transform in the following sense 

j \{ d Ve sin [(2f e + l)ip e ] A(y? e ) < 6 TtVe >, (86) 
J eer 

A(v? e ) < Or, Ve > = I II sin We + W < >> ( 87 ) 



< O r ,r > 

1 Je 



where A(6) = sin 6. 
Proof : - 

Let us prove the first formula. The RHS is written 



RHS = E / U** ( I Il^esin [(2 Je r + l) Ve ] ^S, e (Ge)) II 

{wfff} J e* V eer KG / eir VG 



1 ; 

where G e = Yl e *De9e*- I n the following it is always understood that e e A,e* G A* unless 
otherwise specified. We first write 



sin [(2j e r + l)<p e ] = sm [( 2je + jjjgj x sinv9e = Xj r(^J x A(v? e ), 



sm^ e 



(89) 



and then integrate ^(^e) over V?e using the identity (A8). One obtains 



>,e^r J e* Veer VG J Ve^r G / 

which is 

rhs = w n*, fn^r) n|».(c«) «») 

ie J e* Veer ^ / e VG 

= < O rj r >= LHS, (92) 

and prove the first duality formula. The other one follows from the first one and the orthog- 
onality relation 

oo 

sin(rf^) sm(d j( f)) = ^(5(0 - 0) - 5(6 + 0)). (93) 
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Example : A simple application of this duality formula is given by the case where M = S 3 , 
A being its simplest triangulation with two tetrahedra, one tetrahedra being the interior of 
a 3-ball in S 3 the other one the exterior. We take T to be the 1-skeleton of the triangulation, 
this is a tetrahedral graph. In that case, the gauge fixing PR-observable is the square of the 
67-symbol, 

Jei Je 2 Je 3 



< O r , je > A = f f f , (94) 

^ Je 4 Je 5 Je 6 J 

while the particle PR-observable is related to the Gram determinant 

< Or, Ve >a= 1 • (95) 

A/det(cosy? e ) 

This is the measure for G A /{G x AdG), the duality between these observables was first 
presented in [45]. The above proposition is a generalization of the structure which was 
discovered there. It is also a generalization to the classical case of the duality relation 
proposed by Barrett [46] . 



V. HAMILTONIAN FORMULATION AND TIME EVOLUTION 

In this section we present a hamiltonian formulation of three dimensional gravity which 
generalizes the formalism of Waelbroeck [15] and allows the inclusion of spinning particles 
(more details will be given in [48]). We also briefly discuss its quantization and give an 
Hamiltonian interpretation to the PR amplitude with time. 

A. Classical analysis 

We consider a 2dimensional cauchy surface with punctures and denote A one of its 
triangulation. The graph dual to A is an oriented trivalent graph with open ends ending at 
the punctures We can construct the phase space of 3d gravity by assigning to each oriented 
edge e of T a set of variables (X e , g e ). We call e the internal edges and e v the edges starting 
from an open end v. g e G SU(2) denotes the parallel transport of the spin connection along 
the edge and X e = X£J a is a Lie algebra element 13 characterizing the direction and length of 
the edge of the triangulation dual to e. More precisely, we can obtain X e by the integration 
of the frame field along edge of the triangulation dual to T. X® = j £t dx 11 g* & e a (g*^)~ l where 
g x Se denotes the parallel transport from the starting vertex s e of e to x. The phase space 
algebra is given by 

{X:, g e } = g e J a ; {X e a , X b e } = -2e ab c X c e , (96) 
together with the relations 

g^ e = g- 1 ; X_ e = -g e X e g~\ (97) 
where e denotes any edge and — e denotes the reversed orientation. 



13 J a = J a = icr, where a a are the Pauli matrices, a a <Jb = &ab + i£abc<Jc, is a basis of the Lie algebra 
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Given a face / of T and a vertex sj of / we define the holonomy around / starting from 
Sf, G(f) = I J 9e- We include only internal edges in the definition of G(f), and if / 

J-J-eC/ 

contains an open end v we chose Sf to be the terminal point of e v . We denote G(f) = 
Z7(/)Id + P a (f)J a , note that f/(/) 2 + P a (f)P a (f) = 1 since we are in SU(2). We call 
P(f) the momentum of /. Given an internal edge e C / we can define the 'position' 
Qe(f) and 'orbital momenta' L s (f) of e with respect to / as follows. First, lets introduce 
X g (f) = Gf f (f)Xe{Gs e f {f))^ 1 where G s / f (f) is the parallel transport from Sf to s s along /, 
then 

Q s (f) ee U(f)X g (f)+X s (f)AP(f) + Xsif ^ if) P(f), (98) 
L g (f) ee -(P(/)) 2 X e -( /) + [/(/) (l e -(/) A P(f) + Xg( g (/ [ (/) P(/)j , (99) 

where (X A P) a = e a bc X b P c and X ■ X = X a X a . We can check that 

{01(f), W)} = {LIU), P b (f)} = e\cPV), (100) 

g e -(/)AP(/) = L e -(/). (101) 

Given any internal vertex -u of T we define J(-u) = ^2 e . Se= ^ X e . 

The constraints of the theory are the flatness constraints and the Gauss law 

G(/) = l; J(v)=0, (102) 

for all internal vertices v and all faces which do not contain open ends. These constraints 
are first class. For faces /„ which contains open ends v we impose particle constraints 

9eMfv)9e v = h ev ; X ev = 2S V J . (103) 

where V = 4nm v is half the deficit angle created by the presence of a particle of mass m v , 
S v is the spin of the particle and ho v = exp(9 v Jo)- Among these six constraints, four are 
second class and two are first class, they generate time reparametrisation and U(l) gauge 
symmetry. It is easy to check that P(f v ) and J(v) = —l/2X_ £v + L s (f) commute with these 
constraints. The Dirac bracket involving these observables is therefore equal to the original 
bracket and reads 

{ J», J b (v)} = e ab c J c (v), {J a (v), P b (f v )} = e ab c P c (fv). (104) 

In term of these observables the constraints are 

P a (fv)Pa(fv) = sin 2 6 v ; J a (v)P a (f v ) = S v sine v . (105) 

The Gauss law generates gauge transformations 

{ J A , X e } = [X e , A(s e )]; { Ja, g e } = 9 e A(s e ) - A(t e )g e , (106) 

where A(v) is a collection of Lie algebra elements associated with internal vertices and Ja = 
J a (v)A a (v). The momentum generates translational symmetry. In order to understand 
its action we choose a rooted maximal tree Tr of T passing through all internal vertices, 
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with a distinguished vertex R (the root). We denote by G^(T) the parallel transport along 
the tree from any internal vertex v to the root. We choose a collection of Lie algebra 
element 4>(f) for all faces of T. These elements are all sitting at the root from the gauge 
transformation points of view. We need to transport them to the faces using G§(T). Finally, 
we define P = J2 f pa (f)[( G ? f ( T ))~ 1( f ) (f) G *( T )}a- P<f> commutes with all group elements 
g e , its commutation with X e is given by 

{P^X e } = -J>(/,e) (U(f)$.,{f) + [PM, $.,(/)]) ■ (107) 
where e(f, e) denotes the relative orientations of / and e and 

P Se {f) = (GZifT'PiWZif); kAf) = (G«(T)GZ(f)y\(f)G*(T)GZ(f), 

G(T),G(f) being the parallel transport along T, f. This is the hamiltonian version of the 
translational gauge symmetry defined in (56). We can now introduce a parametric time by 
taking P^ to be our hamiltonian. With respect to this time evolution all J(v)'s X Cv and 
P(/)'s are conserved. However the J(t>)'s or X e ,Q e (f) are not conserved and we could use 
one of them as a clock in order to define a relational time [16, 49]. 

B. Quantization 

The quantization of this system is implemented by taking the spin network representation 
presented in section II C 3 and more precisely in section VII. ®(r,j e ,i v )(9e] 9e v ) £ <8>vVi v 
is a gauge invariant functional of all g e 's valued in the tensor product of representations 
associated with open ends. The dependance on the open end group elements is explicitly 
given by 

$(rj e) j„)(ffel g ev ) = ® v D Iv (g^)<5> [T ^ Iv) (g- e - 1), (108) 

where Dj(g) is the matrix element of g in the representation Vj. 

In the spin network representation the <? e 's are acting by multiplication and the X e 's are 
acting by left invariant derivative 

g e $ { r)(g e ) = 9e^r)(g e ); X:<S> (r) (g e ) = ^ ( r)G? e J a ). (109) 

The physical states are obtained by 'projecting' the spin network states on the kernel of 
the constraints. This 'projector' is rigourously defined as a rigging map [50] mapping spin 
network states onto distributional spin network states (Linear functionals on the space of 
spin networks). This rigging map is easily obtained: We have to multiply the spin network 
states by a factor 5(G(f)) for each face which does not contain open ends. We multiply 
by a factor S(G(f v )g ev hQ v g~^) for each face with open ends. Eventually we have to insert a 
spin projector Pf v j for each open end where Pfj = \S, I) (I, S\, \S, I) being the normalized 
vector of spin S in Vi, (I, S\ being its conjugate. This is clearly a projector, moreover it 
satisfies Pf v v Iv Di v (J ) = —2iS v P T J Iv . By an analysis similar to the one performed in eq. 
(36), one can check that 

P(f v ) a Av)Pf;j v ^(r, je ,i v )(g S ;ge v ) = S v sm9 v P^ Iv <l> {Wv) (g s -,g ev ), (110) 



26 



and leads to the imposition of the spin constraint (105) at the quantum level. The physical 
scalar product between two spin network states is obtained by integrating the product of 
the spin network states with the rigging map over all group elements. It is explicitly written 
in section (VIII). 

One of the main goal is to be able to compute the matrix elements of the unitary operator 
exp Ptf,. We are not there yet, but the Ponzano-Regge time observable we propose is one step 
in this direction. More precisely, lets consider S ffjJl the genus g surface with n punctures. 
Using the gauge fixing we can collapse the triangulation of E Si „ to a triangulation with one 
vertex. 

In this case the dual graph possess only one face and = tr0GV where Gt = 
\Yi=i[ a i-> bi] Ylp=i UphepUp 1 - Let us define the operator 

U\4>\ = j d 9R^-V{iP g -^ 9R )- (HI) 
In the appendix A we prove that 
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Therefore at the level of distributions we have 

tf|*| = W|-d;)*^- (H3) 

The inclusion of the time operator had as an effect to replace the expectation value of S(Gt) 
by Xj{GT)dj ( see the definition of the time operator and also section VIII), therefore it is 
computing the expectation value of the operator 



VI. MASSIVE SPINLESS PARTICLES 

In this section we are going to describe the computation of transition amplitudes for 
massive but non spinning particles. 



A. Definition of the amplitude 

We are considering a manifold M with boundaries 14 DM = Sj TJ E/, A a triangulation 
of M and A* the dual triangulation. We denote by = A* fl E« and Tf = A* n £/ the 
closed trivalent oriented graph obtained by the intersection of the 2-skeleton of A* with 
DM. We label the edges e,, e/ of Tj, Yf by spins ji,jf and consider spin network functional 
^rWsfeJ, $r f j f (ge f ) defined in section (II C 1). 

In our model, we consider that the particles propagate along the edges of the triangulation. 
In order to describe their propagation and interaction we consider "I an oriented graph with 
open ends whose edges are edges of the triangulation A, the open ends of "I all lie in the 
boundary and are identified with vertices of the boundary triangulation A fl Ejj. The 



the boundaries are possibly empty 
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internal vertices of 1 describe the interaction of particles. We label each edge e of 1 by 
an angle 9 e , it is related to the mass of the particle by 8 e = 47rGm e . This correspond to a 
decoration D = (m e , S e ; I Se , It e ,i v ) where S e = I Se = It e = 0. According to our gauge fixing 
procedure we have to choose a maximal tree T* of A* and a maximal tree T of A — 1. 
We define the transition amplitude to be 




This amplitude does not depend on the triangulation A and the choices of trees. It also 
depends only on the diffeomorphism class of 1 [28]. The inclusion of particles amounts to 
insert a particle PR-observable 

6-lfi,Ue,9e-)=l[6 S .,o6o.(Ge). (115) 
eel 

Expanding the delta function and integrating over g e * we can write this amplitude in a 
Ponzano-Regge like form 15 

(^1*^,1 =e(ii^) f n *>) (Ux 3 m)u{^ f; H' 

> VeeT / Ve^m / Veg-1 / t ^ Jet * J ^ J ^ > 

(116) 

where the summation is over all spins belonging to internal edges of A, the spins of the 
boundary edges ej,e/ given by are not summed over and e ti denotes the six edges 
belonging to the tetrahedra t. 

As we have seen in section II C, we can express the transition amplitude in the polarization 
were e or uj is fixed. In order to change polarization we have to add a boundary term 
Sb = — J dM tr(e Aou). The same is possible at the discrete level, we can express the transition 
amplitude as a propagator (jM.ifSe^Je;) depending on the boundary connection: 

G M ^{g- ev g e - } ) = J2 (IHA) *r / j / Ue / )*r i j i (^)(*r /J/ |*r ilJ - i )M ) n, (117) 

(^Tfjf^T^M,! = J dgeAge;G M ^{ge^ge;)^T f ,jj{ge } )^T t , jt {9e t )- (118) 

The term ^rj(s'e) is the discrete analog of exp J dM ti(e A a;). 
We are now going to give justifications for the amplitude 114. 

B. Link with the discretized action 

At the classical level the particle is inserted by adding to the gravity action a particle term 
Sp m = m J dttr(e t Jo). The total action reads St = J d 3 xtr (e(x) A {F(u) + 4GmJo5p(x))), 
where Sp(x) = J dtS^(x — xp(t)). At the discrete level, the curvature term is replaced by 



using recoupling theory and the expansion of Sg s in terms of characters Se e (G) = Xj(hg s )xj(G) (sec 
appendix A) 
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the holonomy G e and the frame field by the Lie algebra element X e , we choose the discrete 
action to be 



where 9 e = 4irGm e . In the limit where the loops around e become infinitesimal and the 
mass is small compare to the Planck mass we recover the continuum action. If we integrate 
over X e we get for each edge of "I a factor S(G e hg e ). This is not yet the term we want, we 
have to remember that when we introduce a particle via the action Sp m we are breaking the 
gauge symmetry at the location of the particle. In order to define G e we have to choose a 
dual vertex st(e) along the face dual to e, the discrete Lorentz gauge group is acting at this 
vertex G e — ► u^ 1 G e u e . If we integrate over the gauge group action we get 



So the insertion of a particle forces G e to be in the same conjugacy class as he e . Note that 
now the momentum of the particle is no longer labelled by a Lie algebra element as in (16), 
but by a group element u e hQ a u~ x . This property is characteristic of a theory exhibiting 
the symmetry of Doubly Special Relativity [43]. The fact that a spinless particle coupled 
to gravity can be inserted by imposing the holonomy around the particle to be in a given 
conjugacy class can be easily seen from our hamiltonian analysis. It was first recognized, at 
the classical level, in [15] and analyzed in great details in [21]. Our amplitudes extend these 
works to the quantum case. 

C. Physical interpretation of the particle observables 

The partition function contains for each edge e of the triangulation a term 5(G e ) which 
ensures that the curvature around e is zero. The insertion of a PR-observable 0^g e removes 
the zero curvature condition and replaces it by the fact that the curvature around e has 
to be in the conjugacy class of hg e . This can be interpreted as the presence of a particle 
of mass m e = 8 e /4irG moving along e, creating a line of conical singularity. For such an 
interpretation to be valid, we have to check that it satisfies a kind of 'mass conservation 
property'. 



Consider a bivalent vertex of the graph. The two edges of this graph carry group elements 
gi and g 2 which are imposed to lie in the conjugacy class labelled by, say, Q\ and 9 2 . In 




(119) 




(120) 




FIG. 2: particle vertex 
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the triangulation this vertex is surrounded by other edges, with group elements g e which 
are imposed to be 1 due to the flatness condition. Now the overall still satisfies the Bianchi 
identity, which in this case reduces on-shell to g\g 2 x = 1 (provided we choose an ingoing 
orientation for the first edge and an outgoing for the second one). This proves that the 
expectation value of the particle PR-observable is zero unless 6\ = 9 2 . If fact, if there is a 
bivalent vertex, this analysis shows that the amplitude (114) will be proportional to Sg 1 (hg 2 ) 
in this case the physical amplitude is the proportionality coefficient. If we naively put equal 
masses at bivalent vertices in (114) we get an infinite result, this is due to the fact that 
there is an extra gauge symmetry associated with the propagation of a particle which is 
time reparametrisation. 

Consider now a trivalent vertex. This time the Bianchi identity around the vertex reduces 

to 

9i = 929s or UxfiQ^ 1 = u 2 h d2 u 2 l u^u^ 1 . (121) 
This can only holds if mi,m 2 and m 3 satisfy the inequality 

cos( ) < cos — < cos( ). (122) 

2 2 2 

This inequality is the analog of the kinematical inequalities of relativistic particles. It is 
shown in the companion paper [28] that this condition actually corresponds to the existence 
of the intertwiner between the associated three representations of DSU(2). 



VII. MASSIVE SPINNING PARTICLES 

In this section we are going to describe the inclusion of quantum massive spinning particles 
coupled to quantum gravity. 

We are considering a manifold M with boundaries 16 DM = Sj A a triangulation 

of M and A* the dual triangulation. We denote by A* = A* n E< and A* = A* n £/ the 
closed trivalent oriented graph obtained by the intersection of the 2-skeleton of A* with DM . 
We denote by Vi (resp. Vf) the vertices of A fl S« (resp. A n £/). Each Vi is at the center 
of a face of A*, for each Vi we choose a vertex Vi lying in the middle of one of the dual 
edges surrounding Vi and draw an edge called e Vi from Vi to Vi, see figure (3). We denote 




FIG. 3: particle vertex 



the boundaries are possibly empty 
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by Ti = A* Uje^jj the union of A* with these edges, this is an oriented graph with open 
ends. We do the same construction for E/ and denote the resulting graph Tf. We color the 
internal edges of Ti by SU{2) spins j Si and the open edges e Vi of Ti by SU(2) spins I Vi . 
The resulting spin network with open ends is denoted (Ti,je. , I Vi ) an d similarly for 

Given such data we construct a spin network functional denoted by $(r*,j g .,i v .) which 
belongs to 

Hr,i v = L 2 - ® v V Iv )/GW) , (123) 

where |e| denotes the number of edges of T,\v\ the number of internal vertices of T, v the 
open ends of P*, and the action of on G' e ' is given by the structure of the graph. The 
functional is explicitly defined as 

$(r Je -,/„)(<7e-; 9e v ) = (®vCv\ ® e - D^(g € ) ® v D h (g Sv )) E <g>„ (V Iv f , (124) 

where Cy are normalized intertwiners for internal vertices. This is a construction similar to 
the one presented in section II C, but adapted to a triangulation. 

These spin network states arise naturally in the hamiltonian quantization of 3D gravity 
with spinning particles as kinematical states satisfying the Gauss law. The operators g e act 
by multiplication and the operators X e as left invariant derivatives, see (109). 

A. Particle graph functional 

The purpose of this section is to describe the construction of a particle graph functional 
which contains the data characterizing the propagation and interaction of particles in the 
Ponzano-Regge model. The data we need are the following: we consider 1 an oriented graph 
with open ends whose edges are edges of the triangulation A, the open ends of 1 all lie in 
the boundary and are identified with the boundary vertices i>j,i>/. 

Each edge e of 1 is labelled by a mass and a spin (m e , S e ); also each starting point of an 
edge e is denoted by s e and labelled by a SU(2) representation V/ Se , each terminal point of 
an edge e is denoted by t e and labelled by a SU(2) representation Vi te , both I Se — S e , I te — S e 
are positive integers and each internal vertex v of 1 is labelled by an intertwiner The 
collection of data (m e , S e ; I Se , I te , i%) is called a decoration and denoted for short D, the 
decorated particle graph is denoted "l^- 

1. Spin projector 

In order to construct the particle graph functional we need to introduce an important 
object: the spin projector IT} r : G — > Vj <S> Vv which plays a key role in the construction of 
the propagator of spinning particles coupled to gravity 

W lv {u) = D\u- x )Pl r D r (u), (125) 

where -D 7 (w) is the matrix element of u in the representation Vj and Pf v = \s,I)(I',s\, 
\s, I) being the normalized vector of spin s in V}, (/, s\ being its conjugate. It is clear that 
I — s,I' — s' E N in order for this definition to make sense. The spin projector is characterized 
by the following properties 

n^OOnj^ii) = 5 s , s ,TL s IJ „(u), (126) 
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n s i r (hu) = Tl^ p^u), for all h £ H, (127) 

n^,(^) = DV^n^H/^'G^for ail «? £ g, (128) 

n?,r(l) = ^> (129) 

The first property justifies the name projector for II, the second property shows that II is 
a functional on H\G. The space H\G is isomorphic to the momenta space of a particle: 
In the Lorentzian case, the momenta p £ M 3 of a 3d particle lie in the upper sheet 7i ~ 
SO(2, l)/SO(l, 1) of the hyperboloid, due to the mass shell condition p 2 = m 2 . Given p we 
can choose a section u : 7i — > G. If a functional f{u{p)) is invariant under the right action 
of H then / doesn't depend on the choice of the section and is therefore a functional of the 
momenta. In our case the space of Euclidean momenta is the sphere, and u is interpreted 
as a particle momenta. The third property shows that II intertwines the action of boosts 
on the momenta and on the external Lorentz index. The last property is a normalization 
condition showing that we are describing a spin s particle. Note that II does not depend on 
the mass of the particle. 

One of the main justification for introducing II is that it appears naturally in the compu- 
tation of the Feynman propagator of spin s particle propagating in flat space. Let s (x) be 
a field of spin s, which is a function valued in V s and let us denote by d^(j) s (x) the traceless 
n-th derivative of s , it is a function valued in V n <g> V s since is specified once we choose 
a traceless symmetric product of vectors in V\. We can project V n ®V s onto V n+S and denote 
by d^ n 4> s \x) the corresponding field. The flat space Feynman propagator for this field is 

(tf'-pW'-'^Hp) = ni ' J,( f p)) . (iso) 

p 2 — m 2 + te 

This is proven in [26] in the three dimensional case, following the analogous proof of [47] in 
the case of dimension 4. 



2. Framed particle graphs 

In order to construct the particle graph functional we need to give a framing to 1. We 
associate to "I £ A the set of edges of A which touche "I but are not in 1 

T n = {e £ Ai | en 1 ^ and e £ 1}, (131) 

where Ai refers to the one-skeleton of A. Each edge of e £ A can be viewed as a dual face /* 
of A*. We define the tube of ~I, denoted to be the set of dual faces corresponding to edges 
in T-|. So e £ T-| iff /* £ T\ The neighborhood of ~I from the point of view of A* is given 
by together with the set of dual faces /*, e £ ~I. is made up of tubes or cylinders 
surrounding the edges of 1 and the topology of is a 2d surface with holes which can be 
viewed as the boundary of a thickening of "I, the holes being at the location of the particle 
(see figure 4). We denote by 1 f and call it a framed particle graph, a graph in A* with the 
following properties: "V is a graph entirely lying in T^, it is topologically equivalent to "I 
and provides ~I with a framing /. The end points of "K are identified with the boundary 
dual vertices Vi,Vf introduced in the previous section. The framing being given by a choice 
of vectors on ~I pointing towards "V". We induce the decoration of ~I on "V and denote the 
corresponding decorated graph by "Tp, see figure (5). For each vertex v of 1 we choose a 
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FIG. 4: tube surrounding a particle worldline 




FIG. 5: The framed graph 



corresponding vertex v* of "V, however not all vertices of "V are of these type there is for 
each edge of "I a vertex st(e) which comes from the intersection of "V with /*, (see figure 
5). Given an edge e of "Id we consider 




FIG. 6: connection between the framed graph and the boundary spin network 

U s lLl Jg:,u e ) = DjJigir^TllAujDjJgl) G Rom(V Ise , V he ) , (132) 

where g* e is the group element corresponding to the holonomy along the dual edges of "V 
going from s* to st(e) and g% corresponds to the holonomy from t* to st(e) . 
We can finally define the particle graph functional to be 

n n / (g e *,u e ) = (® ve iiv\ ®c6Tnf; e)J(e (^*,u e )) G Horn (® Vi Vi v .; ®v } Vi v ^j . (133) 

The PR-observable that corresponds to the inclusion of interacting spinning particles is given 
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by 

H-,/ (j e ,9e*) = Yl^jefi / \\du e 5(G e u e he e u~ l ) ri-,/ (g e * , u e ) . (134) 
eel J V _ / 
This operator is clearly a generalization of the operators we have already considered. For 
instance if we take the decoration such that all S e and all I v are zero we recover the operator 
(115) describing the propagation of a spinless particle without external angular momenta. 
If we take a decoration such that 9 e = 0, we chose I Se = I te and we sum over S e we recover 
the Wilson line operator (84). 

In the case of spinning particles, the spin network states are functional of the 
holonomies g s along internal edges and of the holonomies (<?g„) along the external edges. 
®(ri,h v i Vi )(9e; 9e v ) (resp. ^(r f , jgf ,i Vf )(9e;9e v )) are valued in ® Vi Vi v . (resp ® Vf V Ivf ). We can 
therefore contract these functionals (we denote them $j, $/ for short) with the particle graph 
functional (133) to get a scalar function 

($/ (<?e; 9e v ))|n n /0' e ,^)l^ (9s, 9e J > ■ (135) 

The transition amplitude is given by the following integral 




Note that we integrate only over the internal boundary variables g g and we fix the value of 
g Sv to be unity. This amplitude depends on the choice of the framing. If we make a change 
of framing along the edge e by adding one more twist one sees that the amplitude is modified 
by a multiplicative factor e tev2Sv . 



B. Spin conservation 

It is now easy to check that the expectation value of the PR-observable (134) satisfies 
a mass and a spin conservation property. The argument is similar to the one presented in 
section VI C. Lets consider a bivalent vertex v — t ei — s e2 , this vertex is surrounded by edges 
whose group elements G e are imposed to be unity by the flatness constraint. Therefore, the 
Bianchy identity reduces to 

(91X^91, =« 2 )- 1 Ce 2 ^ 2 , (137) 

where g% ,g% corresponds to the holonomy from v to st(ei), st(e 2 ). G ei is imposed to be 
equal to u~^h-g ei u ei and similarly for G e2 . This implies that h-g ei is conjugated to h^g e2 , so 
9 ei = ±# e2 . This is the mass conservation. If the sign is +, this means that u ei g% = hu e2 g* e 
where h is an arbitrary Cartan group element. The PR observable (134) contains the product 

D he2 {{giX'nhMD^ (g* Se2 )D Itei ((giX^h-M ^ ( 138 ) 

which is equal to 

5s u s 2 D Ite2 ((hu^gtJ-^P^Dj^ (u ei g* Sei ), (139) 
and the spin is conserved. 
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If 61 = —6 2 then u ei g* ei = ehu e2 g*^ where h is an arbitrary Cartan group element and e 



is the Weyl group element 



-1 

1 



When acting on \I,S), Dj{e) gives \I, — S). We see that in this case the expectation value 
of (134) contains a factor t^-Sj and the absolute value of the spin is conserved. This is the 
quantum analog of the remark at the end of section II B, showing that we have to include 
at once both spins S and —S. If we have two incoming particles and one outgoing at an 
interaction vertex the same reasoning implies that |<Si| + IS2I = IS3I 



VIII. COMPUTATION OF TRANSITION AMPLITUDES 

In this section we present explicit computations of transition amplitudes (114,136) in 
the simple cases where we have non interacting spinless and spinning particles on a genus 
g surface. We show that our general definition is free of any infinities in this case and 
reproduces the physical scalar product one expects from the canonical quantization. We 
also compute explicitly the effect of the exchange of two particles and show that the result 
is characterized by the braiding matrix of the quantum group _D K (SU(2)). We also discuss 
the finiteness of our amplitudes for closed manifold. 



A. The physical scalar product and time propagation 

We are now going to specialize to the case of non-interacting particles. We consider a 
simple triangulation of S 0jTl , the sphere minus n points, in terms of 2n triangles and n + 2 
vertices. The triangulation and its dual are drawn in figure 7. The particles are sitting 
inside the small tubes. 
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FIG. 7: Triangulation and dual triangulation of So, n and dual triangulation of So, n x I 



We consider a triangulation A of S 0i „ x I — 1 n , where ~I n is the graph consisting of n 
vertical unbraided segments {x{\ x I. The triangulation we take contains 6n tetrahedra with 
all its vertices in the boundary, the dual triangulation is drawn in figure 7. We now chose 
a gauge fixing condition. Since all triangulation vertices are in the boundary we need to 
gauge fix only the Lorentz symmetry. In figure 8 we crossed the gauge fixed edges and show 
in the RHS the result of the gauge fixing. We can simplify the picture by using the flatness 
condition around dual faces where all edges but one are gauge fixed to eliminate the last 
edge of this face. After this procedure we get a simpler picture drawn in the RHS of fig. 8 
which contains 4n + 2 non gauge fixed edges and 4n + 2 faces. We label as in the figure the 
group elements around a tube p by kp, kp, g Pt i, g p j where the k's label internal edges and we 
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FIG. 8: gauge fixing around a tube: The crossed edges are gauge fixed. 



integrate over them, we also denote by go,i,goj the edges of the boundary going from the 
tube n to the tube 1. The two vertical faces lying on the tube give us two delta functions: 
8(kp (kp)~ 1 )5(kpg p j(kp)~ 1 (gp^i)' 1 ) , the horizontal face closing the tube gives us 5e p (g p j), 
the faces going from the tube p to p + 1, p — 1, • • • , n — 1 give ^(/^(A;^)" 1 ), the face going 
from the tube n to 1 gives <5(£;^goj(^f) _1 #(h )> eventually, there are two faces encircling all 
the tubes one gives <5(go,«<7M • • • <?n,i), the other one gives S(goj). After integration over all 
/c's variables but one, we are left with the propagator 

Gn(g P ,f,g P ,i) = j dk ^^Sikgpjk'^gp^-^Se^gpj)^ S(g lti - ■ ■ g n ,i)S(go,f)S(g ,i)- (140) 

After gauge fixing and taking into account the flatness condition for g , the boundary spin 
network states depend only on (g p )p=i,..., n - The physical scalar product between two spin 
network states is then simply 

/n 
Yldup&fiuphepUp^iiuphepUp^Siuxhe.Ui 1 ■ ■ (141) 
P =i 

This shows that the physical Hilbert space of n particles on a sphere is isomorphic to 
L 2 (SU(2)/f/(l)) ra_1 and a physical state $p(-u l5 •, w n _i) can be viewed as a distributional 
functional $ on the kinematical Hilbert space using a version of the so called 'rigging map' 
[50] 

®p(g P ) = 8(gi ■ ■ -g n ) j ^n^ P 5feM P V( M p) _1 )j $p(«i, •••,Un-i)- (142) 

We recover in this way the scalar product that we obtained in the Hamiltonian quantization 
in section V. Also the physical Hilbert for n particles coupled to gravity is isomorphic to the 
Hilbert space of n particles on the sphere without gravity and can be written as the tensor 
product of n — 1 Poincare representations (S) p ^e p ,o- However, even if the Hilbert spaces are 
isomorphic the representation of operators acting on it is very different, this is exemplified 
in the next section where we compute the action of the braiding operators. 

If we insert the 'time' operator (85) in the computation of the amplitude its effect is to 
replace the 8 function by a character Xj an d we get 

= J l^dUp^fiuphepU-^iiuphepU- 1 )^ XjiuiK^ 1 • ■■u n he n u^)- (143) 
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B. Spinning amplitudes 



We have seen that the amplitude including spinning particles is given by 




We want to specialize to the case where we have n non-interacting spinning particles on the 
sphere using the triangulation described in the previous section. In this case the spin network 
states are functional of the holonomies (g p ) p =i i ... n around the particle and the holonomies 
(u p )p-i > ... n along the external edges. §i(g p ;u p ) (resp. $f(g p ;u p )) are valued in ® v Vj v . (resp 
<8>v f Vi v )• F° r non interacting particles the particle graph functional factorises as a product 
of projectors. It is therefore easy to contract the spin network functionals with the particle 
graph functional (133). The computation of the amplitude is similar to the one made in the 
previous section, if we restrict to the case where $j = $/ it reads 

V U V , Uy) | ®H \lvi S v ))\ 2 . (145) 

This shows that the physical Hilbert space of n particle on a sphere is isomorphic to 
®v'Hs v where 

H s = {Fe L 2 (G)\ F(h e g) = e ms F(g)\Vh e G H}. (146) 

A physical state <3>p(tti, -,u n -i) G <S> v 7~ts v can be viewed as a distributional functional $ G 
<8> v Vi v on the kinematical Hilbert space 

&p(gv',Uv) = 8{9i---9n) (^^duJigyUvho^Uv)- 1 )^ (® V \I V , S v )) $ P (tti, • • • , tt„-i). (147) 
C. Braiding of particles and general surface 

Using these techniques it is possible to compute the expectation value of a more general 
particle graph. If we take M to be the three sphere and consider a closed decorated particle 
graph 1 D . It is shown in [28] that the PR evaluation Z(M, 1c) is equal to a quantum 
group evaluation using D K (SU(2)) a kappa deformation of the Poincare group. We are now 
going to show this explicitly by looking at the exchange two particles. The effect of gravity 
is to give a non trivial statistic to the particles which is governed by a deformation of the 
Poincare group. We will explicitly show this now. 

We consider a triangulation A of So, n x I — "In, where ~l n is the graph consisting of n 
vertical segments, n — 2 of them are unbraided and the segment 1 is crossing the segment 
two so that the particles are exchanged. The triangulation we consider is the same as the 
previous one except around the particles 1 and 2. The corresponding dual complex is drawn 
in the LHS of figure (9). Following the same path as previously we can gauge fixed the 
amplitude and collapse the complex to the LHS of figure (9) where the gauge fixed edges 
are crossed. Only four faces of the collapsed dual triangulation contribute non trivially, its 
evaluation gives 

J c?A;5(^i i/ ^ 2i/ A;" 1 ^ 2i ^i ii )5(A;^ 2i/ A;" 1 ^i ii )5 ei (^ 2i/ )5 e2 (^i i/ ). (148) 
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The effect of the exchange on the scalar product of two spin network is then 

($ f \$i) = J ^Y[dg p Se p (g p )j Qfig^gT 1 , 9i, • • -0n)$i(0i, 92, • • • , g n )S(gi • • • 9n)- (149) 

The effect of the exchange of particles 1, 2 therefore results in the following operator acting 
on the physical Hilbert space 

ft0i,92$p(«l» «2, • • • , «n-l) = $p(MlVMr lM 2, «1, ' ' ' , «n-l)- (150) 

We recognize the action of the R- matrix of _D K (SU(2)) acting on the vector space 7^ o®1~Lg 2 o 
[51]. 

The effect of the exchange of two spinning particles also result in the action of the 1Z 
matrix of D K (SXJ(2)) on the physical Hilbert space 

^(e 1 ,s 1 ),(e 2 ,s 2 )^p(u 1 , u 2 , • • ■ , Un-i) = ^p{uih ei u^ l u 2 , u u ■ ■ ■ , w n _i), (151) 

where P G ® v H Sv - 

In the appendix we present, starting from a triangulation, the computation of the prop- 
agator in the case of a genus g surface. The result is given by the following propagator 

/g 9 
dkSiJlia^b^llSi^ka^k-^S^kb^k- 1 ), (152) 
i=i i=i 

where [a, b] = aba^b^ 1 denotes the group commutator. The general transition amplitude 
for a genus g with n punctures can be deduced and we recover the expected transition 
amplitude. 

/g n g n 

\\dadkX[du p \<$>( a i ,6 i ,Mp/ie p M p 1 )| 2 (5(JJ[ai,6i] JJwp/i^Mp 1 ). (153) 
i=i p=i i=i p=i 

It is interesting to compute the vacua to vacua scalar product 

j v j ' p=i j 

where dj — 2j + 1 is the dimension of the spin j representation. This formula gives the 
volume of the moduli space of flat SU(2) connections [52]. 
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D. Closed manifolds and finiteness of the modified PR model 

In the previous example we have seen that the Ponzano-Regge transition amplitudes are 
explicitly finite after gauge fixing. In the case of closed manifold this is not true in general. 
Even after gauge fixing the Ponzano-Regge partition function is not necessarily finite for 
closed manifold. This is clearly seen if one tries to compute for instance the partition 
function Z{Tj g , n x S 1 ) which can be expressed as the trace of the propagator (153). It should 
be so since such an expression is computing the dimension of the physical Hilbert space 
which is of course infinite. One natural question to ask is whether we can find manifolds for 
which the gauge fixed Ponzano-Regge partition function is finite. One can find necessary 
conditions, for example if a closed manifold M contains a non contractible torus then the 
partition function is infinite even after gauge fixing. 

Lets consider, for instance, a mapping Tori Ej = Ex// ~j where / is a mapping 
class group element ( a non trivial diffeomorphism of E ffi0 ) and the equivalence relation is 
ix {0} = f(x) x {1}. It is then known that if Ey is hyperbolic then the manifold is atoroidal. 
One would expect the gauge fixed Ponzano-Regge partition function to be finite in that case, 
but we haven't checked. Also the evaluation (63) of the partition function shows that the 
invariant is defined by counting the number of flat SU (2) connections. This is very similar 
to the original definition of the Casson- Walker invariant [53] and we therefore expect the 
gauge fixed PR partition function to be finite for integral homology spheres. 

IX. CONCLUSION 

In this paper, we have given a complete treatment of the quantization of 3 dimensional 
Euclidean gravity in the spin foam language including an analysis of local Poincare invari- 
ance, finiteness, topological invariance and the insertion of massive spinning particles. We 
have given a general prescription allowing us to compute quantum transition amplitudes with 
interacting particles. We have introduced the notion of particle graphs functionals which 
generalize the notion of Feynman graphs for theories coupled to three dimensional gravity. 
We have sketched a new hamiltonian treatment of 3d gravity coupled to spinning particles 
and showed that our amplitude prescription computes its physical scalar product. We have 
presented the link between the spin foam quantization and the combinatorial quantization 
of Chern-Simons. 

We feel that our work opens the way to many new develop ements. First, the treatment 
we have presented of the gauge fixing of the symmetries should allow us to tackle the more 
challenging problem of the spin foam quantization of Lorentzian 3d gravity with particles. 
Then, our treatment of particle insertions in spin foam shows a clear link between spin 
foams and Feynman graphs which need to be better understood and eventually generalized 
to higher dimensional gravity. It would also be very interesting to construct explicitly the 
field theory reproducing the amplitude we have given. Our work gives many hints towards 
the answer. The spin foam formalism we have developed here to include particles should be 
naturally extendable to the case where a cosmological constant is present. One also needs to 
develop a better understanding concerning the insertion of time in the quantum amplitude, 
and the corresponding semi-classical interpretation. Eventually, we have to see wether the 
structures we have introduced here to include matter can be exported for the study of 4d 
quantum gravity amplitudes. 



39 



Acknowledgments 

We would like to thank Etera Livine for his help and numerous key discussions. We 
would like to thank Jose Zapata for discussion on the hamiltonian framework, Karim Noui 
and Alejandro Perez for sharing with us some of their results. We are grateful to Bernd 
Schroers, Justin Roberts and John Barrett for the organisation of the ICMS Edinburgh 
workshop (29 June- 5 July 2003) where some of our results were presented. D. L. has 
benefited from a Menrt and Eurodoc grants and the hospitality of PI. 



[1] E. Witten, "(2+l)-Dimensional Gravity As An Exactly Soluble System," Nucl. Phys. B 311, 
46 (1988). 

[2] S. Deser, R. Jackiw and G. 't Hooft, "Three-Dimensional Einstein Gravity: Dynamics Of Flat 

Space," Annals Phys. 152, 220 (1984). 
[3] S. Carlip, "Quantum Gravity In 2+1 Dimensions," Cambridge monographs on Mathematical 

Physics. 

[4] V. Moncrief, "Reduction Of The Einstein Equations In (2+l)-Dimensions To A Hamiltonian 
System Over Teichmuller Space," J. Math. Phys. 30 (1989) 2907. 

S. Carlip, "Notes on the (2+l)-dimensional Wheeler-DeWitt equation," Class. Quant. Grav. 
11, 31 (1994) [arXiv:gr-qc/9309002]. 

S. Carlip, "Six ways to quantize (2+l)-dimensional gravity," arXiv:gr-qc/9305020. 

R. S. Puzio, "The Gauss map and (2+1) gravity," Class. Quant. Grav. 11, 2667 (1994) 

[arXiv:gr-qc/9403012]. 

[5] K. Ezawa, "Classical and quantum evolutions of the de Sitter and the anti-de Sitter universes 
in (2+l)-dimensions," Phys. Rev. D 49, 5211 (1994) [Addendum-ibid. D 50, 2935 (1994)] 
[arXiv:hep-th/9311103]. 

E. Buffenoir and P. Roche, "Two-dimensional lattice gauge theory based on a quantum group," 
Commun. Math. Phys. 170, 669 (1995) [arXiv:hep-th/9405126]. 

A. Y. Alekseev, H. Grosse and V. Schomerus, "Combinatorial quantization of the Hamiltonian 
Chern-Simons theory. 2," Commun. Math. Phys. 174, 561 (1995) [arXiv:hep-th/9408097]. 
E. Buffenoir, K. Noui and P. Roche, "Hamiltonian quantization of Chern-Simons theory with 
SL(2,C) group," Class. Quant. Grav. 19, 4953 (2002) [arXiv:hep-th/0202121]. 
[6] S. Carlip, "Exact Quantum Scattering In (2+l)-Dimensional Gravity," Nucl. Phys. B 324, 
106 (1989). 

[7] C. Meusburger and B. J. Schroers, "The quantisation of Poisson structures arising in Chern- 
Simons theory with gauge group G x g*," arXiv:hep-th/0310218. 

[8] M. P. Reisenberger and C. Rovelli. "Sum-over-surface form of loop quantum gravity", gr- 
qc/9612035, Phys. Rev. D 56 (1997) 3490; 

J. W. Barrett and L. Crane, "Relativistic spin networks and quantum gravity," J. Math. Phys. 
39, 3296 (1998) [arXiv:gr-qc/9709028]. 

J. W. Barrett and L. Crane, "A Lorentzian signature model for quantum general relativity," 
Class. Quant. Grav. 17, 3101 (2000) [arXiv:gr-qc/9904025]. 

R. De Pietri, L. Freidel, K. Krasnov and C. Rovelli, "Barrett-Crane model from a Boulatov- 
Ooguri field theory over a homogeneous space," Nucl. Phys. B 574, 785 (2000) [arXiv:hep- 
th/9907154]. 



40 



M. P. Reisenberger and C. Rovelli. "Spacetime as a Feynman diagram: the connection formu- 
lation", Class. Quant. Grav., 18:121140, 2001; 

A. Perez, "Spin foam models for quantum gravity," Class. Quant. Grav. 20, R43 (2003) 
[arXiv:gr-qc/0301113]. 

[9] L. Freidel and K. Krasnov, "Spin foam models and the classical action principle," Adv. Theor. 
Math. Phys. 2, 1183 (1999) [arXiv:hep-th/9807092]. 
[10] T. Thiemann, "Introduction to modern canonical quantum general relativity," arXiv:gr- 
qc/0110034. 

A. Ashtekar, "Quantum geometry and gravity: Recent advances," arXiv:gr-qc/0112038. 

[11] G.Ponzano and T.Regge, " Semiclassical limit of Racah coefficients" in Spectroscopic and 
group theoretical methods in physics (Bloch ed.), North-Holland, 1968. 

[12] H. Ooguri and N. Sasakura, "Discrete and continuum approaches to three-dimensional quan- 
tum gravity," Mod. Phys. Lett. A 6, 3591 (1991) [arXiv:hep-th/9108006]. 

[13] G. 't Hooft, "Canonical Quantization Of Gravitating Point Particles In (2+l)-Dimensions," 
Class. Quant. Grav. 10, 1653 (1993) [arXiv:gr-qc/9305008]. 

[14] G. 't Hooft, "The Evolution Of Gravitating Point Particles In (2+l)-Dimensions," Class. 
Quant. Grav. 10, 1023 (1993). 

[15] H. Waelbroeck, "Time Dependent Solutions Of (2+1) Gravity," Phys. Rev. Lett. 64 (1990) 
2222. 

[16] H. Waelbroeck and F. Zertuche, "Homotopy Invariants And Time Evolution In (2+1)- 

Dimensional Gravity," Phys. Rev. D 50, 4966 (1994) [arXiv:gr-qc/9401021]. 
[17] H. Waelbroeck and J. A. Zapata, "2 + 1 Covariant Lattice Theory and t'Hooft's Formulation," 

Class. Quant. Grav. 13, 1761 (1996) [arXiv:gr-qc/9601011]. 
[18] L. Cantini, P. Menotti and D. Seminara, "Hamiltonian structure and quantization of 2+1 

dimensional gravity coupled to particles," Class. Quant. Grav. 18, 2253 (2001) [arXiv:hep- 

th/0011070]. 

[19] C. Vaz and L. Witten, "Reconstructing the metric of space-time from Chern-Simons gravity 
with point sources," Mod. Phys. Lett. A 7, 2763 (1992) [Erratum-ibid. A 8, 469 (1993)]. 

[20] E. Buffenoir and K. Noui, "Unfashionable observations about 3 dimensional gravity," arXiv:gr- 
qc/0305079. 

[21] H. J. Matschull and M. Welling, "Quantum mechanics of a point particle in 2+1 dimensional 

gravity," Class. Quant. Grav. 15, 2981 (1998) [arXiv:gr-qc/9708054]. 
[22] H. J. Matschull, "The phase space structure of multi particle models in 2+1 gravity," Class. 

Quant. Grav. 18, 3497 (2001) [arXiv:gr-qc/0103084]. 
[23] P. de Sousa Gerbert and R. Jackiw, "Classical And Quantum Scattering On A Spinning Cone," 

Commun. Math. Phys. 124, 229 (1989). 
[24] P. de Sousa Gerbert, "On Spin And (Quantum) Gravity In (2+l)-Dimensions," Nucl. Phys. 

B 346, 440 (1990). 

[25] F. A. Bais, N. M. Muller and B. J. Schroers, "Quantum group symmetry and particle scattering 
in (2+l)-dimensional quantum gravity," Nucl. Phys. B 640, 3 (2002) [arXiv:hep-th/0205021]. 

[26] L. Freidel, E. Livine and D. Louapre, "Ponzano-Regge model revisited III: The Field Theory 
limit". To appear. 

[27] K. Noui and A. Perez "Three dimensional loop gravity coupled to point particles", to appear. 
[28] L. Freidel and D. Louapre, "Ponzano-Regge model revisited II: Mathematical aspects; relation 

with Chern-Simons theory, DSU(2) quantum group and link invariant. To appear. 
[29] L. Freidel and D. Louapre, "Diffeomorphisms and spin foam models," Nucl. Phys. B 662, 279 



41 



(2003) [arXiv:gr-qc/0212001]. 
[30] Laurent Freidel, " A Ponzano-Regge model of Lorentzian 3-Dimensional gravity ", Nucl. Phys. 

Proc. Suppl. 88 (2000) 237-240 , gr-qc/0102098. 
[31] Stefan Davids, " Semiclassical Limits of Extended Racah Coefficients ", J. Math. Phys. 41 (2000) 

924-943, gr-qc/9807061. 

[32] L. Freidel and D. Louapre, "Ponzano-Regge model revisited IV: Lorentzian 3D Quantum Ge- 
ometry". To appear. 

[33] S. Carlip, "The statistical mechanics of the three-dimensional Euclidean black hole," Phys. 

Rev. D 55, 878 (1997) [arXiv:gr-qc/9606043]. 
[34] H. J. Matschull, "On the relation between 2+1 Einstein gravity and Chern-Simons theory," 

Class. Quant. Grav. 16, 2599 (1999) [arXiv:gr-qc/9903040]. 
[35] L. Freidel and K. Krasnov, "Discrete space-time volume for 3-dimensional BF theory and 

quantum gravity," Class. Quant. Grav. 16, 351 (1999) [arXiv:hep-th/9804185]. 
[36] C. Rovelli and L. Smolin, "Spin networks and quantum gravity," Phys. Rev. D 52, 5743 (1995) 

[arXiv:gr-qc/9505006]. 

[37] A. Ashtekar and J. Lewandowski, "Projective techniques and functional integration for gauge 

theories," J. Math. Phys. 36, 2170 (1995) [arXiv:gr-qc/9411046]. 
[38] J. C. Baez and K. V. Krasnov, "Quantization of diffeomorphism-invariant theories with 

fermions," J. Math. Phys. 39, 1251 (1998) [arXiv:hep-th/9703112]. 
[39] L. Freidel and K. Krasnov, "2D conformal field theories and holography," arXiv:hep- 

th/0205091. 

[40] N. Kawamoto, H. B. Nielsen and N. Sato, "Lattice Chern-Simons gravity via Ponzano-Regge 

model," Nucl. Phys. B 555, 629 (1999) [arXiv:hep-th/9902165]. 
[41] L. Freidel and E. R. Livine, "Spin networks for non-compact groups," J. Math. Phys. 44, 1322 

(2003) [arXiv:hep-th/0205268]. 
[42] J. D. Roberts, "Skein theory and Turaev-Viro invariants," Topology 34 (1995), 771-787. 
[43] L. Freidel, J. Kowalski-Glikman and L. Smolin, "2+1 gravity and doubly special relativity," 

arXiv:hep-th/0307085. 

[44] E. Witten, "Quantum Field Theory And The Jones Polynomial," Commun. Math. Phys. 121, 
351 (1989). 

[45] L. Freidel and D. Louapre, "Asymptotics of 6j and lOj symbols," Class. Quant. Grav. 20, 

1267 (2003) [arXiv:hep-th/0209134]. 
[46] J. W. Barrett, "Geometrical measurements in three-dimensional quantum gravity," Int. J. 

Mod. Phys. A 18S, 97 (2003) [arXiv:gr-qc/0203018]. 
[47] S. Weinberg, "The Quantum Theory Of Fields. Vol. 1: Foundations," Cambridge University 

Press. 

[48] L. Freidel and A. Zapata, "Hamiltonian description of 3D gravity with spinning particles," To 
appear. 

[49] C. Rovelli, "Partial observables," Phys. Rev. D 65, 124013 (2002) [arXiv:gr-qc/0110035]. 
[50] D. Giulini and D. Marolf, "A uniqueness theorem for constraint quantization," Class. Quant. 

Grav. 16, 2489 (1999) [arXiv:gr-qc/9902045]. 
[51] T. H. Koornwinder, F. A. Bais and N. M. Muller, "Tensor product representations of the 

quantum double of a compact group," Commun. Math. Phys. 198, 157 (1998) [arXiv:q- 

alg/9712042]. 

[52] E. Witten, "On Quantum Gauge Theories In Two-Dimensions," Commun. Math. Phys. 141 
(1991) 153. 



42 



[53] K. Walker, "An extension of Casson's invariant", Annals of Mathematics Studies, 126 Prince- 
ton University Press, 1992. 



APPENDIX A: NOTATIONS ON SU(2) 

We use the following notations for SU(2). The elements of the Cartan subgroup H — U(l) 
are represented using 

" e»* 



The Weyl group consists of the identity and the reflection — > h_^. Let us choose a (non 
necessarily normalized) Haar measure on the group, the integration over the group can be 
written as an integration over the conjugacy classes using the Weyl integration formula 

jT dgf(g) = jf ^ Qf ^ fixhex- 1 )^ d6 (A2) 

where H = U(l) denote the cartan subgroup, A(6) = sin 6* and \W\ is the order of the Weyl 
group. 

The representations of SU(2) are labelled by a half-integer j and are realized on the spaces 
yi ~ C 2j+1 . The matrix elements of representations are given by the Wigner functions 
Dii n (g) satisfying the orthogonality property 

_ y 

D 3 mn (g)D^, n ,(g)dg = 8^,—8 m ^8 n ^, (A3) 

Vg being the volume of the group. This relation can be written in terms of convolution 
product for characters 



Xj'(g)Xj(gx)dg = V G 8jji-^—. (A4) 

G a j 

We define the distribution 5<j,(g) to be the distribution forcing g to be in the conjugacy 
class of htj,. It is invariant under conjugation 8^{g) = 8^{xgx~ l ) and normalized by 

/ 8^{g)f{g)dg= / fixh^x'^dx. (A5) 
Jg Jg/h 

We can write this distribution in terms of characters 

S M = T7- Yl Xj(h<i>)Xj(9), (A6) 



V H 

3 

where Vh = 27r is the volume of the Cartan subgroup. The Weyl integration formula imply 
that 

/ d^\<P)8^h B ) = 1. (A7) 
Jh/w 

This means that we can relate this distribution to 8^(9) the delta function on H/W, 
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APPENDIX B: A DELTA FUNCTION IDENTITY 



In this section we prove the following formula 

d 3 xe MXg)^ ± £ ^ = $ n S(±g), (Bl) 



where X = x l <Ji is in the Lie algebra, g = exp(iOn l (Ji) is a SU(2) group element with G [0, ir] 
and n l rii = 1; a 1 are the Pauli matrices, tr<7j<7j = Sij, e(g) = sign(cos6>), and 5(g) is the delta 
function on the group with respect to the normalized Haar measure. First we evaluate 

J d 3 xe MX9) = (2ir) 3 5 {3) (sin 9n). (B2) 
We can use the familiar identity S^(X) = ^rp^p5(|X|) to write this evaluation as 

-^—5(\sme\) = -^—y^5{e-nn)-^—, (B3) 
(sm6») 2 (sm6>) 2 ^ | cos6»| 

2vr 2 



W - 7v2n) + 5(6 - n(2n + 1))) . (B4) 



(sinfl) 2 

The normalized Haar measure on the group is given by 

dg = -d9(sin9) 2 d 2 n, (B5) 

71 

where d 2 n is the normalized measure on S* 2 , therefore we can write 

'<») = 2^5). 5> -**•). < B6 > 

s (-9) = x£wlL 5 ( d -< 2n + 1 ))- ( B7 ) 

Together with (B3) this proves that 

J d 3 xe itTiX9) = 4vr (5(g) + 5(-g)) . (B8) 
A similar computation shows that 

J d 3 xe^e(g) = An (5(g) - 5(-g)) (B9) 

which proves eq(Bl). 



APPENDIX C: EXPLICIT COMPUTATIONS 



In this appendix, we present explicit computations using the gauge fixing procedure we 
introduced. In particular, we compute partition functions for manifolds S s x S* 1 where £ 9 
is the genus g surface, and transition amplitudes and propagator for S 3 x I. We present the 
explicit computations for the genus 1 case and the results for the general genus, which can 
be obtained in the same way. 
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1. Partition function 



First, the results are the following : the gauge fixed partition function for a manifold 
S 9 x S* 1 is given by 

r. 9 9 9 

Z[Z g xS']= / dkdn{a u h) TT dtudbi 5{T] aAa^b; 1 ) TT Sioika^k-^Sibikb^k- 1 ). 

(CI) 

We give the explicit proof in the genus 1 case i.e Si x S* 1 which is the 3 dimensional 
torus. A triangulation is obtained by considering the triangulation of Si with two triangles. 
It induces a decomposition of Si x / into two prisms, each prism can be triangulated with 
three tetrahedra (see figure 10). By identifying the past and future faces, we obtained a 
triangulation of S x x S* 1 with 6 tetrahedra, 12 faces, 7 edges and 1 vertex (see figure 11). 
The 6 tetrahedra (dual vertices) are denoted Pi, I 1: F 1: P 2 , I 2 , F 2 . 





FIG. 10: Triangulation of a prism with 3 tetrahedra. The tetrahedron owning past face is called 
P, the one owning future face is called F, the intermediate one is called I. 



The non-gauge fixed partition function for Si x S* 1 takes the form (44) with 12 group 
elements and 7 5-functions. The 7 dual faces corresponding to the 7 edges are 

F1P1I2P2F2I1, F2P2I1P1F1I2, 11P2P1I2F1F2, P1P2I2F2F1I1, (C2) 
PiP 2 l2F 2 F 1 I 1 , F 2 l\P\l 2 -, PiI\F\Ii, FiF 2 P 2 Pi, (C3) 




The translational symmetry does not need to be gauge fixed since the triangulation 
possess only one vertex. We only have to perform the gauge fixing of the Lorentz symmetry. 
There are 6 dual vertices, we choose to gauge fix to the identity the following 5 group 
elements (see figure 12). 

giiFi, 9hPi, 9hF 2 , 9iiPii 9PiP 2 - (C4) 
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FIG. 12: Gauge fixing on a maximal tree for the Lorentz symmetry on Si x 5 1 : the thick edges 
are gauge fixed to identity 

After this gauge fixing, the partition function can be written as 

Z[Ei x S* 1 ] = J dg FlF , 2 dgp lFl dgp 2F , 2 dg IlF2 dgp 1 i 2 dg Il p 2 dgF 1 i 2 

S{9F 1 Pi gp 1 1 2 9 Pi F 2 gF 2 h)5(g F2 p 2 gp 2 h 9 Pi Pi 9 f ± i 2 )$(gii p> 9 Pi h 9h Pi 9f 2 h ) 

^ {9F 2 ii9Pii2 {9p2Ii9f 1 i 2 )5 (9f 1 f 2 9f 2 P29piF 1 )$ (9f 1 f 2 ) • (C5) 

Solving first the last 8 function gives g FlF2 = 1, then the three ones remaining on last line 
give 

9hP 2 = a = gJ 2Fi , g Plh = b = g F ^ h , g PlFl —k — g^. (C6) 
This leads for the gauged fixed partition function to 

Z[Ei x S 1 ] = [ dadbdk Siaba^b-^Siaka^k'^Sibkb^k- 1 ). (C7) 

Finally we use the measure dfj,(a, b) to take into account the remaining gauge invariance at 
the last dual vertex. 

The previous computation can be systematically generalized to the genus g case. We can 
triangulate S s with 4g — 2 triangles. This induces a triangulation of S 9 x S* 1 with 12g — 6 
tetrahedra and still one vertex. The Lorentz gauge fixing procedure can be systematically 
conducted along the same lines and we obtained the general result announced. 

2. Propagator and transition amplitudes 

The propagator for genus g case is given by 

Gfa, h; a'i, b'i) = j dkS(Y[ aA^X 1 ) JJ S^ka^k-^d&kb^k- 1 ). (C8) 
J i i 

Again we consider only the genus 1 case, i.e the case of the manifold £i x /, where the 
boundaries are a past and a future surface S x . We consider the same way to triangulate 
this manifold than before, with two triangles on each boundary and 6 tetrahedra (except 
that there is no more past/future identification). We first consider the case of the prop- 
agator between boundary connections. Each boundary carries three fixed group elements 
representing the fixed boundary data, and we have 14 group elements living on the internal 
dual edges. There are 4 flatness conditions associated to internal edges of A, and 6 flatness 
conditions around edges of the boundary (see figure 13) 

The dual graph is given by figure 13. We now perform a gauge fixing according to figure 
13. The gauge fixing on the boundary spin-network leave only two group elements. Writing 



46 



FIG. 13: Dual 1-skeleton of the triangulation of Si x I and maximal tree for the gauge fixing 



explicitly the gauge fixed partition function and solving the delta functions, the remaining 
propagator is 

G(a, b; a', b') = J dkdiaba-H-^a'ka^k-^ib'kb^k- 1 ). (C9) 

Recall that we use the measure d/i(a,b) (see (62)) to fix the remaining gauge invariance at 
the last vertex. The gauge fixing procedure we proposed can be generalized in a systematic 
way for the T, g x / case. We obtain for the propagator 

G(a l ,b t ;a' l ,b^}= f dk <5(Q a^V) JJ S^ka^k-^S&kb^k- 1 ). (CIO) 

J i i 
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